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Abstract 

In this thesis we analyze a very simple model of two dimensional quantum gravity 
based on causal dynamical triangulations (CDT). We present an exactly solvable 
model which indicates that it is possible to incorporate spatial topology changes 
in the nonperturbative path integral. It is shown that if the change in spatial 
topology is accompanied by a coupling constant it is possible to evaluate the path 
integral to all orders in the coupling and that the result can be viewed as a hybrid 
between causal and Euclidian dynamical triangulation. 

The second model we describe shows how a classical geometry with constant nega- 
tive curvature emerges naturally from a path integral over noncompact manifolds. 
No initial singularity is present, hence the quantum geometry is naturally compat- 
ible with the Hartle Hawking boundary condition. Furthermore, we demonstrate 
that under certain conditions the quantum fluctuations are small! 
To conclude, we treat the problem of spacetime topology change. Although we 
are not able to completely solve the path integral over all manifolds with arbitrary 
topology, we do obtain results that indicate that such a path integral might be 
consistent, provided suitable causality restrictions are imposed. 



How to read this thesis? 

After completion of this thesis considerable progress has been made and published 
in [1] and [2]. People interested in [1] and [2] can use this thesis as an introduction 
to the subject. The first chapter is written for a very general audience and can be 
skipped by most readers. In reference [1] spatial and spacetime topology change in 
two dimensional causal quantum gravity are formalized in the form of a string field 
theory. Subsequently, in [2] we uncovered a matrix model underlying the string 
field theory. In a forthcoming publication we show why the continuum amplitudes 
can be obtained from a matrix model. One will see that it emerges from a new 
continuum limit of the one matrix model. 
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Introduction to quantum gravity 



1.1 Classical gravity, physics of the large 

Gravity, omnipresent and inescapable... 

Unlike the other three fundamental forces of nature its reach is universal. All 
objects and substances in the universe are sensitive to the gravitational pull. Be- 
sides being mere slaves to the will of gravity, matter and energy also play a more 
proactive role, since everything inside our universe acts as a source for gravity. 

In our daily lives we are only confronted with the passive side of gravity. If we 
jump, the gravitational pull of the earth will inevitably let us fall back down again. 
The only effect by which gravity reveals itself is by dictating the way we move, 
never do we experience our role as sources of gravity. More generic, in no micro- 
scopic or mesoscopic experiment does the gravitational pull between the objects 
play an important role. The reason for this is clear, gravity is an extremely feeble 
force when compared to the other fundamental interactions. Although we usually 
take this fact for granted, it might strike one as strange that the weakest force 
of nature dominates the motion of objects on the scales relevant in our everyday 
life. The fundamental reason behind this is that the source of gravity only comes 
in one flavor, matter and energy arc always positive causing gravity to be always 
attractive. If gravity would have had both positive and negative charges similar 
to electromagnetism, gravity would not have played any role in our everyday lives, 
since it would have been overshadowed by the other forces. 

One of the key insights that enabled Newton to formulate his theory of grav- 
ity is the realization that gravity is not only important at scales familiar from 
our experiences, but it is also the relevant force at solar system scales and even 
beyond. He realized that the motion of a falling apple is similar to the trajectory 
of the moon, both are caused by the tug of gravity. 



1 



2 



1 Introduction to quantum gravity 



Even though Newton's theory was tremendously successful, as it explained the 
motion of the planets with unprecedented accuracy, Einstein felt uneasy. He em- 
barked on a historic quest to construct a more aesthetic description of gravity, an 
endeavor that turned out to be one of the pinnacles of human ingenuity in recent 
modern history. One of his motivations to construct a more elaborate theory was to 
harmonize Newton's theory with the principles of his own theory of special relativ- 
ity Another motivation that greatly influenced Einstein's work when constructing 
his theory of gravity was the equivalence principle, inertial and gravitational mass 
were measured to be the same with remarkable accuracy. From these incentives 
and a few other rather philosophical arguments he developed the theory of general 
relativity. So by combining aesthetic reasoning with known results from experi- 
mental physics he found a geometrical theory that was seen to describe the real 
world. In so doing he extended the validity of gravitational theory to the largest 
scales possible. In particular, general relativity has allowed us to compute correc- 
tions to Newton's theory that are vital for the study of cosmology. Furthermore, 
Einstein's description of gravity has been tested and confirmed to be valid for the 
largest distances, masses and velocities that we can measure. 



1.2 Quantum gravity, physics of the small? 

What about the converse regime? Does gravity really become weaker and weaker 
when we study nature on increasingly small scales? The simple empirical answer 
is yes, the gravitational interaction is so incredibly weak that it has only been 
tested down to millimeter scales. It was found that at these scales Newton's law 
still holds implying that gravity indeed becomes negligible for the extremely small. 
The theoretical expectations are more interesting however. 

We know that the physics of systems on small distances is well described by the 
laws of quantum mechanics. One of the many peculiar features of quantum theory 
is that it connects small and large scales by virtue of Heisenberg's uncertainty 
principle. To probe physics on increasingly small scales one needs progressively 
larger momenta. Since a large momentum implies large energy one expects gravity 
to become very relevant at the very tiny scales, contrary to the naive extrapolation 
of the classical theory. The scale for which the probe gravitational field becomes 
large is the Planck scale Lp ~ 1.6 x W^^m. At this scale the energy needed to 
resolve the microstructure needs to be so concentrated that a black hole would 
form. 

Quantum mechanics is, unlike gravity, a theory of probabilities. We know that 
all other forces and all matter fields satisfy its probabilistic laws, so why should 
gravity be an exception? To avoid the coexistence of classical and quantum theory, 
gravity should be quantized too. 
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Why have we not yet been able to accomphsh this? What sets gravity apart 
from the other forces of nature? What makes it so hard to unify gravity with 
the laws of quantum mechanics? The reasons are plentiful, one essential fact that 
makes the analysis of general relativity very hard in general, also on the classical 
level, is that it is highly nonlinear. This nonlinearity is much more severe than 
in the other interactions of the standard model as the action is not even poly- 
nomial in its fundamental field, the metric. This has dramatic consequences for 
the quantization of gravity by perturbative methods. For the most straight for- 
ward methods it basically implies that there is an infinite number of interaction 
vertices that have to be taken into account ^. Therefore wc can conclude that 
the standard formulations of gravity are not easily treated by perturbation theory. 
Another glaring problem that exemplifies the tension between gravity and pertur- 
bation theory is the absence of a natural dimensionless coupling constant to define 
a perturbative expansion. Instead, the coupling constant of gravity, Newton's con- 
stant Gn, has dimensions of inverse energy squared. Subsequently, the natural 
parameter of the perturbation expansion is G^E'^ and we see that the coupling of 
gravitons increases with their energy. At the point where the gravitons reach the 
Planck energy the coupling becomes strong, GnE"^ ~ 1, which inevitable leads to 
a breakdown of the perturbation expansion. In contemporary terms we say that 
the dimensionful nature of Newton's constant makes general relativity nonrenor- 
malizable as a quantum field theory. 

Several points of view can be taken regarding the nonrenormalizability of gravity. 
The most popular stance is that the problem comes from an inherent mismatch 
between the principles of general relativity and quantum theory. According to 
this attitude, a resolution for a quantum theory of gravity can only be found in a 
modification of the physical principles behind cither quantum mechanics, general 
relativity or both. The most popular candidate for such a scenario is string theory, 
where gravity is found to be compatible with quantum theory only if it is accom- 
panied by a plethora of extra fields and dimensions. Gravity by itself is viewed 
as a mere low energy effective theory, and the exact harmonization with quantum 
theory happens only upon considering the dynamics of the fundamental strings. 

A perhaps more conservative attitude is to suppose that gravity and quantum 
mechanics are not fundamentally incompatible per se, but that the standard per- 
turbation theory simply is an inadequate tool for the quantization of gravity. Pre- 
cisely this philosophy is an inspiration for the models we present in the present 
thesis. 



^With a suitable (non tensorial) field redefinition of the metric it is however possible to write 
the Einstein action in a polynomial form. This implies that it is possible do perturbative quantum 
gravity with a finite number of interaction vertices. 
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The construction of a nonperturbative formulation of gravity is a far from trivial 
task however. For example, even though the other field theories of the standard 
model are considerably simpler than gravity, we cannot solve the path integrals 
exactly. Often, the path integrals are merely a helpful tool to set up a perturba- 
tive description of the physical problem at hand. Although extremely successful 
in QED, it is not an adequate scheme to study physics in strong coupling regimes 
such as confinement in QCD. The computation of quantities that go beyond per- 
turbation theory is often very difficult, but even worse, it is mostly unclear whether 
there exists a nonperturbative definition of a path integral at all! In more than two 
dimensions there are few methods that enable one to address the nonperturbative 
existence of path integrals. In most cases the nonperturbative definition of a path 
integral in field theory is only possible by defining it as a limit of a discrete theory. 
Although mathematically more rigorous, it is in practice not a very convenient 
tool to compute concrete amplitudes, analytical methods are largely unavailable. 
Nonetheless, the enormous growth in computing power over the recent years has 
transformed lattice quantum field theory from a mathematically nice idea into a 
serious competitor in the arena of theoretical physics. In particular, the study of 
QCD has benefitted a lot from these developments. Among the successes are the 
calculation of realistic values for meson and baryon masses from first principles. 
Such formidable achievements are currently beyond reach of other methods. 

In this thesis we investigate simple gravitational models that are based on the 
method known as Causal Dynamical Triangulations. In spirit the scheme is a suc- 
cinct gravitational analogue of lattice QCD. It is a natural method to define the 
path integral by a lattice regularization. What remains is a finite statistical sum 
that, similar to lattice QCD, lends itself perfectly to computer simulations. One 
distinguishing feature that sets Causal Dynamical Triangulations apart from other 
discrete attempts is that a genuine causal structure is imposed on the quantum 
geometry from the outset. The results of the simulations are encouraging, in four 
dimensions a well behaved continuum limit seems to exist and there is compelling 
evidence that a classical spacetime superimposed with small quantum fiuctuations 
emerges from the nonperturbative path integral. Despite the intriguing results 
these numerical methods have to offer, the understanding is far from complete 
and inherently restricted by computer power. Furthermore, the statistical model 
is very complicated and has so far resisted attempts at a solution by analytical 
methods. 

In two dimensions the situation is much better however, the pure gravity model 
can be explicitly solved and many interesting results can be obtained. Of course 
one might contest that two dimensional gravity is an oversimplified model as it 
does not possess some of the essential difficulties of four dimensional gravity such 
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as a dimensionful coupling constant. Nevertheless it still contains some vital char- 
acteristics that set gravitational theories apart from any other. Issues such as 
diffeomorphism invariance, background independence and the Wick rotation are 
as relevant for the two dimensional model as they are for its higher dimensional 
analogues. 

Besides being interesting from a pure quantum gravity point of view, it might 
also be regarded as a minimal form of string theory. Particularly, the two dimen- 
sional model of Causal Dynamical Triangulations might shed some light on the role 
of causality on the worldsheet of a string. A tantalizing indication that this might 
indeed be consequential is that the results of two dimensional Causal Dynamical 
Triangulations are physically inequivalent to the outcomes of two dimensional Eu- 
clidean quantum gravity. 

For the purposes of this thesis we primarily view two dimensional quantum gravity 
as an interesting laboratory were nonperturbative aspects of quantum gravity can 
be studied in an exactly solvable setting. Before presenting our original contribu- 
tions, wc first discuss some general remarks and present the known results of two 
dimensional Causal Dynamical Triangulations in chapter 2. 

In chapter 3 we start the discussion of our first generalization by reviewing the 
previously established relation between Euclidean and Causal dynamical triangu- 
lations. It is discussed that imposing causality has the important consequence that 
the spatial topology of the geometries in the path integral is fixed. Additionally, 
we recall that in Euclidean Dynamical Triangulations, as for example defined by 
matrix models, the quantum geometry is highly degenerate in the sense that the 
spatial topology fluctuations dominate the path integral. 

In the remaining sections of chapter 3 we show that this situation is not as black 
and white as is discussed above. In an original contribution we demonstrate that 
one can allow for spatial topology change in two dimensional causal quantum grav- 
ity in a controlled manner. Wc argue that the topology fluctuations are naturally 
accompanied by a coupling constant reminiscent of the string coupling. Upon 
taking a suitable scaling limit we show that the quantum geometry is no longer 
swamped by the topology fluctuations. Surprisingly, we are able to compute the 
relevant amplitudes to all orders in the coupling and sum the power series uniquely 
to obtain an exact nonperturbative result! 

In chapter 4 we return to the "pure" model of Causal Dynamical triangulations. 
In this chapter we extend the existing formalism by studying boundary conditions 
that lead to a path integral over noncompact manifolds. We begin by recalling 
that a similar mechanism is familiar from non-critical string theory where the non- 
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compact quant.Tim geometries are known as "ZZ brancs". Further we show that 
a space of constant negative curvature emerges from the background independent 
sum over noncompact spacetimes. Fascinatingly, we can compute the quantum 
fluctuations and are able to show that they are small almost everywhere on the 
geometry! The model is a nice example of how a classical background can appear 
from a background independent theory of quantum gravity. 

To conclude, we tackle the problem of spacetime topology change in chapter 5. 
Although we are not able to completely solve the path integral over all manifolds 
with arbitrary topology, we do obtain some results indicating that such a path 
integral might be consistent, provided suitable causality restrictions are imposed. 
As a first step we generalize the standard amplitudes of causal dynamical triangu- 
lations by a perturbative computation of amplitudes that include manifolds up to 
genus two. Furthermore, a toy model is presented where we make the approxima- 
tion that the holes in the manifold are infinitesimally small. This simplification 
allows us to perform an explicit sum over all genera and analyze the continuum 
limit exactly. Remarkably, the presence of the infinitesimal wormholes leads to 
a decrease in the effective cosmological constant, reminiscent of the siippression 
mechanism considered by Coleman and others in the four dimensional Euclidean 
path integral. 
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2D Causal Dynamical 
Triangulations 

As explained in chapter 1 there is as yet no satisfactory theory of four dimensional 
quantum gravity, even though both quantum mechanics and general relativity have 
been formulated over eight decades ago! Many obstructions to the unification of 
the two theories, both technical and conceptual, still remain after all this time. 

2.1 Quantum gravity for D < A 

Part of the complications of quantum gravity disappear in lower dimensional mod- 
els for quantum gravity, making them interesting playgrounds where one is not 
confronted with all the issues at once. So one can consider the study of lower 
dimensional models as an action plan to tackle the problems step by step. Of 
course such a simplification comes at a price, some vital features of the real-world 
four dimensional theory are lost. The salient property that distinguishes the four 
dimensional theory from its lower dimensional analogues is that it possesses two 
propagating degrees of freedom whilst the lower dimensional theories have none, 
a fact that can be shown by a canonical analysis. Despite missing this essential 
characteristic, there is a host of problems that lower dimensional models still share 
with the four dimensional theory. An example is the dimensionful nature of the 
gravitational coupling constant in both three and four dimensional gravity. Con- 
sequently, both theories are perturbatively nonrenormalizable by power counting. 
From a canonical analysis however, it is known that there are no local degrees of 
freedom so one deduces that there are at most finitely many degrees of freedom. 
For a comprehensive review of three dimensional quantum gravity see [3]. If four 
dimensional quantum gravity shares the same dissimilarity between the pertur- 
bative and nonperturbative descriptions, it might also be a much better behaved 
theory than the perturbative expansion leads us to believe. 

Often the impression is created that since three dimensional quantum gravity only 
contains finitely many degrees of freedom by canonical analysis, the theory can 
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be completely solved. This is a deceptive representation of the state of affairs 
though. There is no single model that is generally accepted by the theoretical 
physics community. Many problems are unsolved and different approaches give 
different results for important conceptual problems such as, do space and time 
come in discrete units or not? Another issue that does not yet have a satisfactory 
explanation within three dimensional quantum gravity is the explanation of black 
hole entropy of the BTZ black hole [4]. An interesting proposal to do this was 
recently put forward by E. Witten [5] where he basically defines the gravity theory 
by its two dimensional boundary conformal field theory and relates the black hole 
entropy to the degeneracy of states of that conformal field theory. The article is 
also a nice example of the fact that three dimensional quantum gravity has not 
yet been solved in all details and that points of view keep changing as Witten 
also personally changed his viewpoint on the subject. In a seminal work [6] he 
showed that the theory could be written as a Chern-Simons gauge theory which is 
a fairly simple gauge theory that can be quantized and he was of the opinion that 
this equivalence should also hold in the quantum theory. Now on the contrary, he 
advocates that the equivalence is only valid semiclassically since, amongst other 
issues, the Chcrn-Simons formulation docs not require the viclbcin to be invertible 
whereas the metric formulation does. His current opinion is that Chern-Simons 
theory is a useful tool, only to be used for perturbative arguments and it is not rich 
enough to fully capture all aspects of three dimensional gravity such as the physics 
of black holes. Of course a lot of these statements rest on opinions and conjectures 
and, although interesting, one should treat them with care. The argumentation 
with respect to black holes for example rests on the assertion that three dimensional 
black holes are a pure gravity phenomenon. This statement can be called into 
doubt for the obvious reason that black holes usually form by collapsing matter 
distributions. The situation in higher dimensional gravity is more complicated 
though since it seems to be possible to form a black hole by collapsing gravitational 
radiation. Subsequently, three dimensional black hole physics might not tell us 
something about pure gravity but it might describe gravity coupled to matter. So 
we conclude that three dimensional quantum gravity is an interesting arena where 
a lot more can be said than for four dimensional gravity. Many of the fundamental 
issues that it shares with the four dimensional theory remain unsolved however. 
The next step down the ladder is two dimensional quantum gravity. In this step we 
lose part of the perturbative analogy to four dimensional theory, in two dimensions 
Newton's constant is dimensionless and the Einstein Hilbert action becomes a 
topological term, making the theory renormalizable by power counting. Even 
though it is even less similar to four dimensional gravity than the three dimensional 
theory in this respect, it still possesses important conceptual characteristics such as 
background independence, diffeomorphism invariance and the problem of defining 
a Lorentzian theory. A very appealing advantage of working in two dimensions is 
that there exits a plethora of exactly solvable models, both discrete combinatorial 
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as continuum, that can be treated with techniques from conformal field theory and 

statistical mechanics. 

Beyond being a toy model for four dimensional quantum gravity, the two dimen- 
sional model is also interesting from the string theory point of view. To discuss 
this relation let me give a sketch of some of the basic principles behind string the- 
ory. The most convenient way to define string theory is to start from the Polyakov 
action. It essentially describes the string as two dimensional quantum gravity cou- 
pled to scalar fields [7], where the scalar fields act as the embedding coordinates 
of the string. In string theory however, the two dimensional world sheet metric 
is introduced as a mere auxiliary variable. If one uses the equations of motion 
the Polyakov action reduces to the Nambu Goto action which is written purely 
in terms of the coordinates and metric of the embedding, or equivalently, the tar- 
get space. The reason why the Polyakov action is at least locally equivalent to 
the Nambu-Goto action is that it is invariant under Weyl rescaling of the world 
sheet metric. The world sheet metric in general has three independent compo- 
nents, where two can be seen to be gauge degrees of freedom from diffeomorphism 
invariance and the third is unphysical because the action is invariant under Weyl 
transformations. 

Therefore, classically the Polyakov and Nambu-Goto formulations are equivalent 
but quantum mechanically this is not true in general since the measure of the path 
integral over world sheet metrics is not invariant under Weyl transformations. 
Commonly this property of the quantum theory is referred to as the conformal 
anomaly [7, 8, 9]. Subsequently, in general the conformal factor of the world 
sheet metric is not a pure gauge degree of freedom implying that the quantum 
theory based on the Polyakov action is not locally equivalent to a, so far unknown, 
quantum theory employing the Nambu-Goto action. One can however remedy this 
situation by coupling precisely 26 scalar fields to the world sheet in which case the 
conformal anomaly is precisely cancelled. This suggests that the bosonic string 
naturally "lives" in a 26 dimensional target space^. 

To gain a deeper understanding of string theory, people also investigate the dy- 
namics of strings in dimensions different from 26 where the conformal anomaly 
is not cancelled by the target space scalar fields and it has to be interpreted in 
a different way. The study of these models is appropriately dubbed non-critical 
string theory. In this language pure two dimensional quantum gravity is referred 
to as c = non-critical string theory (see for example [10]), where c is a quantity 
called 'the central charge' and is related to the expectation value of the trace of the 
energy momentum tensor of the matter fields coupled to two dimensional quantum 
gravity. 

"^As is well known, bosonic string theory is unstable and possesses a tachyon. To resolve this 
problem one needs to add fermions and supersymmetry. In the resulting superstring theory the 
target space manifold is 10 dimensional. Upon considering nonperturbative effects it is expected 
however that the theory should be described by membranes embedded in a 11 dimensional target 
space. 
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In this thesis we focus on exactly solvable two dimensional gravity models. We 
mainly regard these two dimensional gravity models as a testing ground for higher 
dimensional models for quantum gravity. Nevertheless, since spatial sections of 
two dimensional spacetimes are in fact one dimensional objects we switch between 
(non-critical) string theory and gravity terminology depending on the application. 

2.2 Problems and solutions in quantum gravity 

Taking two dimensional gravity seriously as a model for quantum gravity one has 
to deal with some of the same issues that one faces in the quantization of four 
dimensional gravity. Let me highlight some of the fundamental questions that one 
faces in any background independent approach: 

1. How should one deal with the problem of time in general relativity? Can 
one define a notion of time and define a Wick rotation? 

2. As in any gauge theory one is instructed to factor out the volume of the 
gauge group to avoid divergencies. So in the case of gravity one is faced with 
factoring out the group of diffeomorphisms. Can we do this in any practical 
way and is it possible to find a regularization procedure compatible with this 
symmetry group? 

3. Which class of geometries should be included in the path integral? Should 
the spatial topology be fixed? Should one include geometries with arbitrary 
spacetime topology? 

One of the oldest and most influential ideas to deal with question 1 is due to Hawk- 
ing who takes the pragmatic point of view that one should start with a Euclidean 
formulation from the beginning [11]. The hope was that once the Euclidean theory 
was solved one would be able to find a natural Wick rotation. Of course ignor- 
ing problem 1 simplifies the quantization procedure to some extent but still four 
dimensional Euclidean quantum gravity shares many of the problems such as 2 
and 3 with its Lorentzian counterpart. Hence the hope of solving four dimensional 
Euclidean quantum gravity and performing a Wick rotation afterwards has so far 
not materialized into a concrete theory. 

The quest to answer question 2 has been somewhat more successful. In [12] Regge 
realized that if one introduces a specific lattice regularization one can formulate 
the dynamics of classical general relativity without explicitly referring to a partic- 
ular coordinate system [13]. An intensely studied model in four dimensions that 
utilizes Regge's ideas is quantum Regge calculus. In this approach the topology 
of the lattice is fixed and the length of the edges are the fundamental dynamical 
degrees of freedom [14]. Although an interesting approach it has met with some 
technical difficulties that have kept the theory from providing clear cut results on 
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the nonperturbative sector of the quantum theory. Particularly, it is not clear how 
to define the measure, several proposals exist but there does not seem to be a 
general consensus. 

To retain the benefits of Regge's coordinate invariant geometry but at the same 
time avoiding some of the technical issues associated with quantum Regge calculus, 
the method of dynamical triangulations was developed (see [15] for a comprehen- 
sive review). In this approach the same lattice regularization is used as in Regge 
calculus, but instead of fixing the lattice and promoting the edge lengths to dy- 
namical variables, the lengths of the edges are fixed and the lattice itself becomes 
the dynamical object. Unlike Regge calculus, dynamical triangulation methods 
are not optimally suited to regularize a given smooth classical geometry but are 
highly efficient methods for defining a measure for the path integral over geome- 
tries. Dynamical triangulations are particularly effective in two dimensions where 
the models reduce to systems that can be exactly solved by methods known from 
statistical physics. Especially the use of matrix models and their large N limit 
[16] turned out to be particularly fruitful for the construction of two dimensional 
Euclidean quantum gravity models, see for example [17, 18]. 

It has been shown that the results of these dynamical triangulation models coincide 
nicely with results from continuum calculations in the conformal gauge as intro- 
duced by Polyakov [7]. He showed that the dynamics of two dimensional gravity 
can be obtained from a nonlocal action, often called the induced action. In the 
conformal gauge this action is equivalent to the Liouville action, which is a local 
action. About a decade ago the interest in the field of two dimensional Euclidean 
gravity was revived since it was shown in two seminal works [19, 20] that Liouville 
theory can be quantized using conformal bootstrap methods. 
Note that the induced action does not represent any local propagating degrees 
of freedom, in accordance with canonical considerations. Consequently, Liouville 
theory also does not describe any local metric degrees of freedom either, since 
it is a gauge fixed version of the induced action. Liouville theory does however 
provide nontrivial relations between global geometric properties of the quantum 
geometries such as the volume and the length of its boundaries. We would like to 
stress that it implies that two dimensional gravity is not topological, at least not 
in the strict sense of the word since it depends on the metric information of the 
manifolds. The term topological is not used unambiguously however, one example 
is Witten's Chern Simons representation of three dimensional gravity. Often this 
theory is referred to as a topological theory while in fact it does encode metric 
information in an explicit fashion. The reason for this confusion is twofold. Firstly, 
the Chern Simons theory does not encode any nontrivial local metric degrees of 
freedom but only describes global characteristics of the manifold related to the 
metric. Secondly, one does not need a metric to write Chern Simons actions in 
general, so a Chern Simons theory where the gauge field does not represent any 
metric degrees of freedom is a topological theory. A second example where the 
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word topological is not used in the strictest sense is topological string theory. Here 
the term topological is used to indicate that the dynamics of a string is insensitive 
to the local geometry of the target space but as in gravitational Chern Simons 
theory the global metric properties of the manifold are important. 
Because of the exact solvabihty of the two dimensional Euchdean models, they 
realize the first step in Hawking's attitude to quantum gravity in the sense that 
they are explicit solutions of Euclidean path integrals. Despite the analytical 
control one has however not been able to take the next step. An unambiguous 
continuation to Lorentzian signature has so far not been found. 
The successes of dynamical triangulation methods for two dimensional Euclidean 
quantum gravity inspired J. Ambj0rn and R. Loll to develop a dynamical triangu- 
lation theory that also addresses question 1, known as Causal Dynamical Triangu- 
lations (CDT) [21]. The idea behind the CDT approach is that the path integral 
should only contain histories that have a built in causal structure. The suggestion 
that one should enforce causality on individual geometries in the path integral 
goes back at least to Teitelboim [22, 23]. In CDT the triangulations are given a 
definite causal structure by imposing a particular time slicing and a fixed spatial 
topology. A fundamental distinction with respect to the Euclidean models is that 
in CDT one considers discretizations of spacetimes with a genuine Lorentzian sig- 
nature. Given the time slicing one can make a clear distinction between timelike 
and spacelike edges which allows one to define a Wick rotation that converts the 
quantum mechanical sum over probability amplitudes into a weighted statistical 
mechanical sum. The statistical model that one obtains after applying the Wick 
rotation has been exactly solved for the two dimensional model and encouraging 
results have been obtained for three and four dimensional models using computer 
simulations [24, 25, 26, 27, 28, 29]. 

As discussed above, the results from the Euclidean dynamical triangulation models 
are corroborated by continuum conformal gauge calculations. Similarly, the results 
from the two dimensional CDT model can also be obtained by a continuum calcu- 
lation. Even before the advent of CDT, Nakayama showed that in the proper time 
gauge, two dimensional quantum gravity reduces to a simple quantum mechanical 
model [30]. Using this fact he derived the same amplitudes that one obtains in 
CDT. So if the Euclidean model is equivalent to quantum gravity in the confor- 
mal gauge and CDT is related to quantum gravity in the proper time gauge, one 
would expect that the results of the two theories coincide, since they merely reflect 
a different choice of gauge. How can it be that calculations in the two different 
gauges leads to different results? Is gauge invariance broken? On the continuum 
level these questions are not completely understood, but on the dynamical trian- 
gulations side this problem can been analyzed in detail [31]. From this analysis it 
is clear that the Euclidean path integral contains many more geometries than the 
CDT. One intuitively sees that the Euclidean path integral contains geometries 
where the proper time gauge cannot be chosen globally. 
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Question 3 on the issue of topology change is a highly debated and controversial 
topic, the possible answer seems to vary immensely from approach to approach. 
In most conservative approaches to quantum gravity the stance is taken that one 
should first figure out the quantization of gravity on a manifold of fixed topology 
and only a posteriori consider the possibility of topology change. Even though 
this statement seems rather unambiguous it creates a bifurcation between meth- 
ods that are inherently Euclidean by nature and methods that take Lorentzian 
aspects of gravity seriously. In most methods incorporating some Lorentzian as- 
pects one makes the additional assumption that also the spatial topology is fixed. 
For theories based on Euclidean geometry there is no a priori distinction between 
space and time, implying that a fixed topology of just space might not be very 
natural. 

In more radical theories such as Group Field Theory (GFT) [32] the point of view is 
very different since the change of topology of space and time are an essential ingre- 
dient in its formulation. An even more radical view is taken in for example causal 
set theory, a theory where causality is elevated to the main guiding principle [33]. 
In this approach the concept of a manifold is abandoned from the beginning and 
replaced with points that possess an elementary causal ordering. Consequently, 
the configuration space of causal sets is exclusively contains by topological rela- 
tions. 

Recapitulating, the two dimensional CDT model is one of only very few exactly 
solvable models known to the author that addresses both questions 1 and 2. One 
of the foremost objectives of this thesis is to build on the success of this strategy 
and to present models where we address all three of the questions posed above. In 
particular, we take the process of spatial topology change into account in this ex- 
plicitly Lorentzian setting by introducing a coupling constant for this interaction. 
Luckily we can make a detailed analysis of the sum over spatial topologies, since 
we are able to solve the model to all orders in the coupling constant and sum the 
series uniquely to obtain a full nonperturbative result for this process! In chapter 
5 we also address the issue of spacetime topology change from two different angles. 
Although we are not able to obtain the same level of nonperturbative control as 
for spatial topology changes, interesting results are obtained. 
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Before going into the details of the dynamical triangulation approach we discuss 

a particular continuum aspect of the path integral that we wish to compute. In 
particular, we discuss the role of boundaries in a gravitational path integral. 
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2.3.1 Boundaries and preferred frames 

Although all explicit path integral calculations in this work are performed in the 
context of 1 + 1 dimensional quantum gravity we start our discussion in the setting 
of 3 + 1 dimensional gravity. By presenting the arguments regarding the time 
variable in four dimensions we emphasize that the issues are as relevant for real 
world 3 + 1 dimensional gravity as for our two dimensional models. 
Formally, we define the path integral as, 

Ga{X,Y;T) = J V[g]e-'^^\ (2.1) 

where the S[g] is the standard action for general relativity for manifolds with 
boundaries. So the approach we take to quantum gravity is a rather minimal 
one. No extra degrees of freedom beyond the metric are introduced as is done for 
example in string theory and we do not need a non standard action containing a 
new undetermined parameter such as the Barbero-Imirzi parameter that appears 
in the Hoist action [34]. This action is the classical starting point of present day 
loop quantum gravity. The most familiar form for the action of general relativity 
for manifolds with boundaries was introduced by Gibbons and Hawking [35] and 
reads as follows. 



S[g:Gr,] = -—- / (fx^gR+—— / d^xVhK. (2.2) 

The first term is the standard Einstein-Hilbert action and the second term is the 
Gibbons-Hawking- York boundary term. The boundary term is introduced to make 
the variational principle well defined for manifolds with boundaries. In effect, this 
term cancels the second derivatives in the Einstein Hilbert action such that one 
does not need to specify the derivatives of the metric, only the metric itself needs to 
be given. This fact is particularly clear in the first order, or equivalently Palatini, 
formalism where we can explicitly write the Gibbons-Hawking- York boundary term 
as a total derivative, 

S{e,uj,GN) = —^ [ eabcd{e''Ae'AF'''-d{e''Ae'Au;'"')). (2.3) 

Here F'^'^ is the curvature of local Lorentz transformations and is defined in terms 
of the spin connection, 

= rfo;"^ + a;", A a;^^ (2.4) 
Combining the bulk and boundary contributions in one term gives 

S{e, uj, On) = / Cabcd (e^A e^A u\ A u^" + 2e» A de^ A u^) . (2.5) 

IbTTGjv Jm 
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Note that this action does not depend on derivatives of the spin connection! In 
the first order formahsm the vielbein and the spin connection are treated as inde- 
pendent fields. Without couphng to matter one obtains the following equation of 
motion for the spin connection, 

eabcd (e"A e''A 00% + e" A de^S^) = 0. (2.6) 
Using this equation the action can be written concisely as, 

S{e, u, Gm) = T-rV / '-^-^ ^ ^ ^'') ■ (2.7) 

leTTGAT Jm 

Contracting equation (2.6) with an epsilon symbol gives the more familiar form of 
the equation of motion of the spin connection, 

eabcde^AT^^O, (2.8) 

where T'^ is the torsion two form, 

= {de^ + w\ A e^) . (2.9) 

So the equation of motion of the spin connection yields the constraint that the 
connection is torsion free. Together with the tetrad postulate, 

d,e\ + uj;,e\-Vl,e\^Q, (2.10) 

this fully determines the affine connection and the spin connection in terms of 
the vielbein. The affine connection can be expressed in terms of the metric only 
and reduces to the familiar Levi-Civita connection. By again invoking the tetrad 
postulate one can derive the following form for the spin connection 



^J■ 



2^ e u;^«;A 



= e"'^e^" {n^^x - ^nx^ + ^^Amk) • (2.11) 



Here ^ij,kx is the object of anholonomy 

^i^kX = e;i5[,e,A], (2.12) 

which one can view as objects that measure how much the tetrad basis deviates 
from a coordinate basis. Combining (2.7) and (2.11) we can obtain the action 
defined entirely in terms of the vielbein, 

S{e,GN) = -r^ ! e,,ed(e"Ade''Aa;^(e)). (2.13) 



Written in this form it is clear that the action only depends on first derivatives 
of the vielbein which is what we wanted to show. It should be stressed that this 
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is merely a rewritten form of the conventional Einstein-Hilbert-Gibbons-Hawking- 
York action, no new physics is introduced. 

An important point about the gravitational action including boundary term (2.13), 
is that it is no longer invariant under local Lorentz transformations! This is evi- 
dent, since the spin connection does not transform tensorially under local Lorentz 
transformations. So including boundaries in a gravitational path integral neces- 
sarily introduces a preferred Lorentz frame. The existence of a preferred Lorentz 
frame should however not be confused with concepts as diffeomorphism invariance 
and background independence. These notions are not incompatible with the exis- 
tence of a preferred Lorentz frame. Furthermore, one should realize that the choice 
of a specific Lorentz frame does not affect the bulk dynamics as the local Lorentz 
transformations are still gauge transformations for the bulk geometry, only the 
boundaries break the symmetry. 

2.3.2 Summary 

In the previous subsection (2.3.1) we recalled the role of boundaries in classical 
general relativity. It was stressed that the action for a manifold with boundaries 
is not invariant under local Lorentz transformations. In other words, boundaries 
naturally introduce a preferred Lorentz frame. Within the discrete framework of 
causal dynamical triangulations the preferred frame of the boundary is used to 
define a time foliation of the manifolds in the path integral. 

2.4 Topology change 

In this thesis we discuss some models where we lift the constraint on the (spatial) 
topology and allow for geometries that have handles and/or baby universes in a 
constrained way. More concretely, we present models where a coupling constant 
is introduced for the splitting of a string. If one inchides these more complicated 
geometries in the path integral, the spatial sections of the geometry have the 
topology of several S^^s. This means that one is in fact considering a multi-particle, 
or better multi-string, theory. History has shown that the framework of quantum 
field theory is the best way to deal with multi-particle quantum theories. So in 
our case the best way to deal with the baby universes and the handles would be 
to develop a (non-critical) string-field theory. Although a full-fiedged string-field 
theory based on CDT is beyond the scope of this thesis, we do show that even 
nonperturbative results in the coupling constant for the string interaction vertex 
can be obtained! 

In section 3.2.1 we argue that manifolds with spatial topology change do not admit 
a Lorentzian metric everywhere. However, in the models we discuss the Lorentzian 
signature of the metric only vanishes at a countable number of points. The analysis 
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of the Wick rotation around such points we leave to future work, for the present 
purposes we confine the discussion of our models to the Euclidean domain. 
In chapter 3 we present results where we perform the sum over all tree diagrams 
of our interacting non-critical string theory based on CDT. In particular we obtain 
the disc function and propagator that are nonperturbatively dressed with string 
interactions. In chapter 5 we go beyond tree level and investigate the loop expan- 
sion, enlarging the class of geometries to include manifolds of arbitrary spacetime 
topology i.e. arbitrary genus. The genus expansion is however considerably more 
complicated than its tree level counterpart, hindering us to find nonperturbative 
results in the coupling. So our analysis is limited to perturbation theory, which is 
used to compute results up to order two in the genus expansion. 
Furthermore, in chapter 5 we make an attempt to go beyond perturbation theory 
even when considering the genus expansion. Even though we are unable to sum 
the genus expansion in all generality we can study some non-perturbative effects in 
a toy model where we constrain the holes to stay at the cutoff scale. Limiting our 
focus to quantum geometries that satisfy this constraint one can allow the number 
of holes to be arbitrary and we can perform the sum over genera explicitly. An 
interesting implication of the model is the suppression of the vahie for the effective 
cosmological constant, reminiscent of the suppression mechanism considered by 
Coleman and others in the context of the four dimensional Euclidean path integral. 
The study of higher genus manifolds and baby universes in the path integral ap- 
proach to quantum gravity has received considerable attention in the context of 
two dimensional Euclidean quantum gravity. As has been mentioned, two dimen- 
sional Euclidean quantum gravity is a different quantum theory of 2D gravity 
where the metrics in the path integral have Euclidean signature from the outset. 
Note that if one applies the Wick rotation as described above to the causal prop- 
agator it is also defined by a path integral over Euclidean geometries. The set of 
Euclidean geometries in the causal propagator is however only a very small subset 
of the geometries included in the path integral for Euclidean quantum gravity. 
This can be understood by observing that in Euclidean quantum gravity one does 
not enforce the topology of each spatial universe to be an S^. The relation be- 
tween causal and Euclidean quantum gravity can be made precise if one defines 
their respective path integrals by dynamical triangulation methods. It turns out 
that the Euclidean theory is precisely related to the causal theory by removing 
baby universes [31]. One can also show the reverse relation by starting with causal 
dynamical triangulations and then adding baby universes [21]. 
The fact that the Euclidean theory contains baby universes and the Causal theory 
does not, leads to the conclusion that the Euclidean theory is strictly speaking 
not a theory of one single string whereas two dimensional causal quantum gravity 
is, if one views both as non-critical string theories. Although this is an appealing 
way to view the dynamics of the string one must be careful since it is a picture 
that is purely based on the worldsheet geometry, the relation to the dynamics 
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in target space is not a priori clear. Interestingly, there is no weight associated 
with the branching of baby universes in the Euclidean theory which allows them 
to proliferate and actually dominate the path integral in the continuum limit. 
This leads to the pecuhar situation that the dynamics of non-critical string theory 
defined through Euclidean quantum gravity is largely independent of the dynamics 
of an individual string but is dominated by the multi-particle, or "multi-string", 
nature of the theory^. Specifically, it can be shown by direct calculation that 
at each point of the quantum geometry there is an outgrowth, or equivalently a 
baby universe, at the scale of the cutoff. One of the prime consequences of this 
dominance of baby universes is the non canonical dimension of time, exemplifying 
the fractal nature of its quantum geometry. 

The model presented in chapter 3 could be viewed as a theory where the quantum 
geometry is allowed to form baby universes arbitrarily as in Euclidean quantum 
gravity with the essential difference that we introduce a coupling constant for each 
baby universe. It is shown that this weight effectively tames the proliferation of 
baby universes preventing the amplitudes to be dominated by cutoff scale out- 
growths. The predominant signal that illustrates the mechanism is the canonical 
scaling dimension of time, which is intimately related to the fact that the Haus- 
dorff dimension is two as in the case of "pure" CDT [21, 36] and not four as is the 
case in the Euclidean theory [37, 38, 39]. 

2.5 Simplicial geometry 

In this section we describe how to nonperturbatively define and compute path 
integrals in two dimensional quantum gravity by the method of causal dynamical 
triangulations. One of the principles on which the method is based, is the fact 
that most quantum field theories can only be defined beyond perturbation theory 
by implementing a lattice regularization. 

2.5.1 Quantum particle from simplicial extrinsic geometry 

A classic example where the lattice regularization plays an important role in the 
definition of the path integral is the non relativistic free particle. Although a 

very simple system, it shares some of its essential features with path integrals for 
quantum gravity. Before going into the details of the gravitational path integral 
we first discuss the quantization of the non relativistic particle in some detail and 
highlight the features that are similar to the gravitational quantum theory. One of 
these similarities is that both systems are examples of "random geometry" i.e. the 
individual histories in the path integral have a geometrical interpretation. A slight 
difference however is that in quantum gravity the histories only contain information 
about their intrinsic geometry whilst the dynamics of the non relativistic particle 



^for more accurate assessment of this statement see (3.1) 
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is determined by the extrinsic geometry encoded in its velocity as the classical 
action is given by 

S[x{t)] = J dtx^. (2.14) 

Recall that the central object in the path integral formulation of quantum me- 
chanics is the propagator or Feynman kernel G{x,t;x',f). It describes the time 
evolution of wavefunctions in quantum mechanics 

i;{x,t) = J dx'dt'G{x,t;x',t')ip{x',t'). (2.15) 

Physically, it gives the probability amplitude to measure the particle ai x,t given 
the initial location x',t'. In principle the propagator can be found by deriving it 
from the Schrodinger equation. Feynman showed however that one can instead 
compute the propagator by computing a weighted integral over all continuous 
paths between the initial and final point. 

J P[x(i)]e^^[^W]. (2.16) 

Note that (2.16) merely is a formal expression, one needs to make sense of what 
it means to integrate over all trajectories. To explicitly define the path integral 
for the free non relativistic particle one is instructed to first discretize the space of 
paths to reduce the path integral to a finite dimensional integral. Commonly this 
is done by decomposing a general trajectory of the particle into N piecewise linear 
segments that correspond to infinitesimal time intervals e = {t" — t')/N (fig. 2.1). 
The evaluation of the path integral now amounts to computing N integrals of the 
following form, 

N-l „ f N-1 ^ 

G{x, t; x\ t') = lim A'^ / dxk explij^ ^i^j+i - ^j) \ ■ (2-17) 

fe=i I j=o J 

Note that this integral is not well defined, since the integrand is a complex valued 
phase factor. To be able to compute the integral one has to analytically continue 
the time variable t ^ r = it. This so-called Wick rotation converts the path 
integral into a set of real gaussian integrals that can be performed to obtain the 
regularized Euclidean amplitude. The continuum amplitude can now be computed 
by taking the limit where the time intervals become infinitesimally small, e — > 
0. Notice that in this limit the velocity and therefore the extrinsic geometry 
becomes singular at each point of a typical trajectory in the path integral. In 
other words, the typical paths that contribute to the path integral are highly non- 
differentiable. Similarly, the intrinsic geometry will typically also be singular if one 
defines a gravitational path integral by dynamical triangulations. After taking the 
continuum limit one can obtain the physical Lorentzian amplitude by applying the 
inverse Wick rotation t —>■ t = —it. 
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Figure 2.1: Illustration of the path integral for a one-dimensional non-relativistic 
quantum mechanical problem, e.g. a propagating particle. One possible path of the 
configuration space (path space) is drawn. The "virtual" particle is propagating 
from xq to xn in a piecewise linear path of N steps of time e = {t" — t')/N each. 



Notice that in this form the Wick rotation is an analytic continuation in the 
coordinate t implying that this procedure is not invariant under coordinate trans- 
formations. Taking a gravitational viewpoint, one might ask the question whether 
it is possible to define a Wick rotation for point particles that does not involve 
coordinates? In the following we suggest that it is indeed possible to construct a 
Wick rotation for a point particle that is independent of the particular background 
and coordinate system. 



2.5.2 The invariant Wick rotation for particles 

We start by considering the standard action principle for the massive relativistic 
point particle. The action is well known and simply proportional to the invariant 
length of the four dimensional wordline, 

S = -m j ds, (2.18) 

where the invariant length is given as usual by 

ds^ = -g^^{x)dx^dxr (2.19) 

Suppose that a classical trajectory is written as a;^(r), where r is an arbitrary 
coordinate that labels the points along the worldline. Then the action may be 
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written in the familiar form 




(2.20) 



This action has the very important property that it is invariant under both four 
dimensional coordinate transformations and reparameterizations of the coordinate 
on the worldlinc. Thus it really describes the extrinsic geometry of the worldline 
as it is embedded in the ambient space and not some particular choice of coordi- 
nates. The action has its shortcomings however, since it has a complicated square 
root dependence and is not valid for massless particles. These difficulties can be 
overcome by introducing a one dimensional einbein on the worldline. In terms of 
the einbein the action for a point particle can be written in the following form 

(e"'*' - m^) e dr. (2.21) 

Solving the equation of motion for the einbein gives 

+ eV = 0. (2.22) 

Substituting this result back in (2.21) one recovers the action (2.20). The here 
described relation between the two point particle actions (2.20) and (2.21) is com- 
pletely analogous to the relation between the Nambu Goto and the Polyakov string 
actions respectively. One way to view the einbein is that it is just the lapse func- 
tion of the one dimensional geometry along the worldline, since the einbein has 
just a single component. Essentially, the action (2.21) describes the physics of a 
point particle as one dimensional quantum gravity coupled to the coordinates and 
metric of the target space. Now we are in a position to define the coordinate and 
background independent Wick rotation for the point particle. If we introduce the 
signature constant ^ for both the lapse function of the target space and the lapse 
function of the worldline we obtain the following action. 




(2.23) 



where 

i;2 ^ ^nH"^ + 2S^ix" + KpX^x". (2.24) 

From (2.23) and (2.24) it is now easily seen that the Lorentzian action is rotated 
to i times the Euclidean action when ^ is analytically continued from —1 to 



2.5.3 Quantum gravity from simplicial intrinsic geometry 

As motivated in section (2.5.1), regularization by lattice methods can be a powerful 
tool to define and evaluate geometrical path integrals. In the example of the 
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Figure 2.2: A patch of flat space represented by a regular triangulation 



point particle the path integral is evaluated with the help of a discretization of its 
extrinsic geometry. Similar to the point particle, quantum gravity is a quantum 
theory of geometry, yet unlike the quantum particle it is the intrinsic geometry that 
plays a central role. So the strategy of causal dynamical triangulations to formulate 
a path integral for quantum gravity is to employ a suitable discretization for the 
intrinsic properties of the space of geometries. 

Discrete methods have had a long tradition in geometry and are a natural tool in 
the study of the gravitational dynamics. Notably, Regge [12] discussed that the 
intrinsic geometry of a manifold can be discretized by piecewise flat geometries. 
In arbitrary dimensions a general piecewise flat geometry consists of flat building 
blocks called polytopes. These polytopes are the natural generalizations of polygons 
(d=2) and polyhedra (d=3). Often the piecewise flat geometries are constructed 
solely from elementary polytopes known as simplices. A simplex is the higher 
dimensional generalization of a triangle, in any dimension it is the polytope with 
the minimal number of boundary components. In principle a discrete manifold that 
is constructed purely from simplices is called a simplicial manifold. Often however 
such geometries are simply referred to as triangulations. These triangulations can 
be thought of as a straightforward analogue of the piecewise linear trajectories 
that appear in the construction of the path integral for the point particle. 
One of the beneflts of using a simplicial geometry to approximate an arbitrary 
manifold is that the metric properties are completely fixed by specifying all edge 
lengths. From the metric information one can extract the local curvature which 
is the relevant object when considering gravitational dynamics. The notion of 
curvature for simplicial manifolds was found initially by Regge [12] and was later 
refined in [40]. The prime focus of this work is on two dimensional quantum 
gravity, so we explain the Regge curvature in this simple setting. 
A simple way to understand the Regge curvature is to first consider a regular 
triangulation of equilateral triangles (fig. 2.2). Such a triangulation is a proper 
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simplicial representation of an everywhere flat manifold. This can be seen by 
noticing that each vertex is associated with precisely 6 triangles. Furthermore, the 
triangles are equilateral implying that the angle between its sides is ^n. So wc can 
conclude that the total angle around each point is equal to 27r as should be for a 
flat manifold. Introducing local curvature deformations can now be done in two 
different ways. Historically, the preferred approach is to deform the triangles by 
altering the length of the edges. Applying such a deformation to one of the edges 
of the flat triangulation causes the total angle around its vertices to be different 
from 271. Since all triangles are still constrained to be flat this implies that conical 
singularities are introduced by the deformation. The amount by which the total 
angle around a vertex differs from 2n is called the deflcit angle e^, = 2n — ^j-,^ 6*^. 
The scalar curvature of a vertex can now be directly related to this deficit angle 
by invoking the concept of the so-called dual lattice, namely 




where K is the volume of a cell of the dual lattice. Each triangulation has a 
unique dual lattice that is easily constructed by connecting the "baryccnters" of 
the triangles with edges of the dual lattice (flg. 2.3). The relation between the 
scalar curvature and the deflcit angle is deflned by parallel transporting a vector 
along the edges of the dual lattice that encompass the vertex under consideration. 
The curvature is proportional to the angle by which the vector is rotated after 
one full encircling of the cell of the dual lattice. The result is given by (2.25), the 
angle by which a vector is rotated after one revolution along the dual lattice is 
proportional to the strength of the conical singularity and inversely proportional 
to the volume of the cells of the dual lattice. 

So the curvature of a geometry which is discretized according to the prescription 
of Regge is given by a set of Dirac delta functions located at the vertices of the 
triangulation. The distributional character of the curvature is very similar to the 
non differential behavior of the histories in the path integral of the point particle 
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Figure 2.4: A positive deficit angle created by deforming a triangle 

and related to the distributional properties of quantum fields in general. 

Given the above notion of curvature it is straightforward to introduce the discrete 

equivalent of the Einstein Hilbert action, called the Regge action, 



Observe that this discrete action is manifestly independent of any coordinates as 
is promised by the title of [12]. Originally this discrete formulation of the gravi- 
tational action was conceived as a useful tool to study classical aspects of general 
relativity. It was proposed that the gravitational dynamics could be conveniently 
studied by varying the length of the edges of the simplicial manifold. This ap- 
proach to simplicial gravity is called Regge calculus. In its classical incarnation it 
has had some success but the interest particularly gained impetus with the propo- 
sition that it might serve as a convenient platform for constructing a quantum 
theory of gravity [14] [41]. The main idea is to use the Regge action to construct 
a path integral over the edge lengths of a simplicial manifold with fixed connectiv- 
ity. The approach is succinctly called Quantum Regge calculus, since it is a rather 
faithful generalization of Regge calculus to the quantum domain. 
Quantum Regge calculus has not been able give us much understanding beyond 
semiclassical gravity however. In the context of two dimensional gravity some 
doubts have been raised about the consistency of the approach [42] . One of the 
objections is that the approach is not able to reproduce results from other methods 
such as Liouville field theory or dynamical triangulations. 

In [43] it was proposed that although the formalism of Regge calculus is free of 
coordinates it is not completely gauge invariant. According to their point of view 
it is possible to perform such a gauge fixing, but the associated Faddeev-Popov 
determinants generate a highly non-local measure which makes the theory very 
hard to handle. 

An insightful way to visualize the possible overcounting problems in Regge calculus, 
at least in two dimensions, is to consider the discretization of a flat two dimensional 
manifold in terms of four squares (fig. 2.5). If the squares are equilateral it is 
clear that the total angle around the central vertex is 27r, hence the manifold 




(2.26) 
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Figure 2.5: All discretizations for which + x + /5 + 7 = 27r represent the same 
flat geometry. 



is flat everywhere also at the central vertex. If we allow the edges connecting 
the vertex to fluctuate there still are some possibilities that keep the manifold 
flat. Basically, the fluctuations of the edges that keep the total angle around 
the central vertex fixed at 27r might be considered "gauge transformations" , since 
they do not alter the intrinsic geometry of the manifold. In this context it is 
interesting to notice that the Regge action only depends on the total angle around 
a vertex and not on the angles of individual simplices. One might however argue, 
that the overcounting problems only appear when one discretizes manifolds with 
a high degree of symmetry, such as flat space. Moreover, in the path integral 
these special geometries are typically of "measure zero" which implies that the 
overcounting issue is only a minor problem. 



2.5.4 Dynamical triangulations 

To ameliorate the technical issues of the Regge calculus program, the method of 
dynamical triangulations was developed. Similar to Regge calculus the scheme is 
free of coordinates, since both methods are based on the Regge action. The crucial 
difference however, is that in dynamical triangulations the length of all edges in 
the triangulations are fixed and the geometry is encoded in the nontrivial gluing 
of the simplices. So instead of altering the edge lengths one introduces curvature 
by adding or removing simplices (see fig. 2.6). An indication that this method 
does not lead to overcounting problems in simple situations is that the simplicial 
representation of a given flat manifold is unique. In addition, adding or removing 
a triangle always changes the physical curvature and volume of the geometry, 
implying that also in more complicated situations the discretization seems to be 
free of overcounting problems. Although dynamical triangulation methods are 
not very efficient to approximate individual classical geometries, they are ideally 
suited for a discretization of geometries as histories in a path integral for quantum 
gravity. 
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Figure 2.6: A positive deficit angle created by removing a triangle 



The dynamical triangulation approach to quantum gravity was initially introduced 
in the context of two dimensional Euclidean quantum gravity [44, 45, 46] where it 
turned out to be a very powerful technique to explicitly compute the continuum 
amplitudes for Euclidean quantum gravity. Contrary to the usual situation, the 
discrete approach of dynamical triangulations is often more powerful than contin- 
uum methods. Many amplitudes can be obtained with ease and the procedure is 
considerably more straightforward than the computation of the amplitudes from 
the corresponding continuum Liouville field theory. Actually, the quantization 
of Liouville theory only gained considerable impetus about 15 years later in the 
seminal work [19]. Although the results do not yet cover all aspects of quantum 
Liouville theory, it is an important contribution to the quantization of two dimen- 
sional gravity. Another example of the strength of the dynamical triangulation 
method is that the results for the amplitudes for Euclidean manifolds of arbitrary 
genus are exclusive to this method [47, 48]. 

So the approach of dynamical triangulations is proven to be a useful tool to study 
Euclidean quantum gravity. It turns out however that if one studies the quantum 
geometry of two dimensional quantum gravity it does not behave as we might 
expect from a "realistic" quantum theory of gravity. Not even the dimension of 
geometry is what it is supposed to be, the Hausdorff dimension of the quantum 
geometry is four and not two. The principle cause of this behavior comes from the 
domination of cutoff scale outgrowths in the path integral. In chapter 3 we analyze 
the properties of the quantum geometry of two dimensional Euclidean quantum 
gravity in some detail. 

Since the results of the dynamical triangulation scheme can be compared with 
continuum calculations, one concludes that the fractal structure of the geometry 
is not a consequence of the triangulation method, but an intrinsic property of two 
dimensional Euclidean quantum gravity. The initial attitude was that the some- 
what degenerate behavior of the quantum geometry is due to the simplicity of 
two dimensional gravity. Therefore, higher dimensional versions of the dynamical 
triangulation model were developed and investigated by means of computer simu- 
lations, see [49] and [50] for the three dimensional model and [51, 52] for the four 
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dimensional model. 

The hope that the higher dimensional models might be better behaved than the 
two dimensional model did not materialize though. It was found that in the infinite 
volume limit these higher dimensional dynamical triangulation models have two 
phases, a crumpled phase where the Hausdorff dimension is very large and a tree 
like phase where the geometry resembles a branched polymer. Both phases are 
not satisfactory from a physical point of view, so a suitable continuum limit is not 
automatically reached in either phase. Nonetheless, one of the initial ideas was 
that an appropriate continuum limit perhaps exists at the critical point separating 
the crumpled and branched polymer phases. Subsequent analysis revealed those 
hopes to be in vain as it was shown that even in the four dimensional model the 
phase transition is of first order [53, 54]. 

An important lesson that can be learned from these investigations is the follow- 
ing, if one constructs a random geometry model based on simplices of a certain 
dimension one is not at all assured that the quantum geometry behaves anything 
like a manifold of that dimension. 

2.6 2D causal dynamical triangulations 

In this section we introduce the concept of causal dynamical triangulations. The 
motivation behind the inception of causal dynamical triangulations was twofold: 

1. The quantum geometry of four dimensional Euclidean dynamical triangula- 
tion models seem incompatible with a well behaved continuum limit. 

2. The Lorentzian signature of the spacetimes should be taken seriously. 

To address the Lorentzian nature of the path integral it is of paramount importance 
to have an intrinsically defined notion of time. In causal dynamical triangulations 
this problem is addressed by studying piecewise linear geometries that have a 
layered structure. The layered structure of the triangulations allows one to globally 
distinguish timelike and spacelike edges. Furthermore, the global foliation of the 
discrete geometries allows one to define a consistent Wick rotation. 
The central amplitude one aims to compute in two dimensional causal dynami- 
cal triangulations is the so-called cylinder amplitude or causal propagator. This 
quantity describes the probability amplitude for a quantum ensemble of two di- 
mensional geometries of Lorentzian signature with an initial and a final boundary, 
where every point of the initial boundary has the same timelike geodesic distance to 
the final boundary. In two dimensions the boundaries are one dimensional curves 
with the topology of a so the only information that characterizes the intrinsic 
geometry of the boundaries is their length. 
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Figure 2.7: Section of a 2d Lorentzian triangulation consisting of spacetime strips 
of height At = 1. Each spatial shce is periodically identified, such that the simplicial 
manifold has topology [0, IJxS"^. One sees that a single strip with lower boundary 
length It and upper boundary length It+i consists exactly of It up pointing triangles 
and It+i down pointing triangles. 



The strategy of causal dynamical triangulations is to discretize the manifolds in 
the path integral with fiat triangles. These Minkowski triangles have one spacelike 
edge satisfying = and two timelike edges with = —a^. By construction, 
one chooses spacetimes which consist of T/a strips with the topology of x [0, 1], 
where the timelike height of the strip is proportional to a. Each strip has two 
spatial boundaries, one spatial section at time r with length L{t) = alt and one at 
time T + a with length L{t + a) = alt+i- The geometry of a general circular strip 
is now determined by the ordering of It+i triangles "pointing up" and It triangles 
"pointing down" as illustrated in fig. 2.7. 

Since the triangles are genuine patches of fiat Minkowski space, they are naturally 
equipped with a local light cone structure. Furthermore, from the global distinction 
between timelike and spacelike edges one sees that the triangulation equips the 
manifold with a global causal structure. By virtue of this global causal structure 
it is possible to define a Wick rotation for curved manifolds. In the triangulation 
context the Wick rotation amounts to changing the squared length of the timelike 
edges from negative to positive signature = —a^ ^ Li^ = a?. As in the 
continuum case this rotation should be treated with some care and one must show 
that the Lorentzian action defined with = —a? and the Euclidean action with 
Li = 0? can be connected by a smooth deformation. In appendix A we discuss 
the Regge action for these two dimensional Minkowski triangulations. It is shown 
that the Lorentzian and the Euclidean Regge action, multiplied with the imaginary 
unit i, are indeed connected by a smooth analytical continuation of the parameter 
a from —1 to 1 where a is defined by = aa?. So the Wick rotation indeed 
possesses the desired property that the weight in the path integral is converted 
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from a complex phase factor to a real Boltzmann type weight. 

yy. giSRegge{r'<"-) ^ e"^^'^**''^^'"). (2.27) 

As we shall see in the forthcoming section the above defined kinematical struc- 
tures reduce the computation of the path integral for the causal propagator to a 
statistical mechanics problem. 



The basic ingredients of any statistical model are the entropy and the Boltzmann 
weight. For the two dimensional causal dynamical triangulations model the en- 
tropy is generated by the number of geometrically distinct ways one can organize 
the Minkowski triangles in a layered structure as depicted in fig. 2.7. The Boltz- 
mann factor on the other hand is related to the Rcgge action associated to the 
triangles. In two dimensions the Regge action is particularly simple, since the two 
dimensional Einstein- Hilbert term is a topological invariant. 



where x(M) = 2—20 — 6 is the Euler characteristic of the manifold M, q is the genus 
and b is the number of boundary components of the manifold. For the moment 
we fix the topology of the manifold to be of the form x [0,1]. This choice 
corresponds to what we refer to as "bare" or equivalently "pure" causal dynamical 
triangulation model as originally introduced by Ambj0rn and Loll in [21]. One 
of the new contributions of this thesis is that we go beyond the assumption of 
fixed topology and in chapter 3 we introduce a model where controlled spatial 
topology changes are an integral part of the quantum geometry. For such models 
the Einstein-Hilbert part of the action is essential as it tames the topology changes 
by introducing a Boltzmann weight that suppresses manifolds with complicated 
topology in the path integral. Since the topology is fixed in the pure model the 
Einstein-Hilbert action can act at most as an overall phase factor in the path 
integral. Furthermore, for the pure causal dynamical triangulations model we 
are solely interested in manifolds with the topology of a cyhnder so the Euler 
characteristic is zero, making Newton's constant irrelevant for two dimensional 
quantum gravity without topology change. 

The only quantity which is relevant for the dynamics of the pure dynamical tri- 
angulation model is the total volume of the manifold. The Regge action is then 
simply proportional to the added volume of all triangles of a specific configuration. 
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(2.28) 



^Regge(T) = Aa2Ar(T) 



(2.29) 



where A is the bare cosmological constant and N{T) the number of triangles in the 
triangulation T. Note that an order one factor coming from the volume term has 
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been absorbed into A. The path integral for the propagator can now be written as 
follows 




(2.30) 



where T denotes the causal triangulations with initial boundary length li and final 
boundary length I2, and C-r denotes the volume of the automorphism group of a 
triangulation. Basically it is the symmetry factor of the manifold that is still left 
after factoring out the diffeomorphisms. After a Wick rotation the discrete sum 
over quantum amplitudes is converted to a genuine statistical model with a real 
Boltzmann weight, 

G,{hM,t) = J2^e-''^'^^^\ (2.31) 

where it should be noted that A and A differ by an order one constant because 
of the different volume of Minkowskian and Euchdean triangles (see appendix A) . 
The layered structure of the triangulations has the natural implication that the 
propagator satisfies the following semi-group property or composition law, 

Gxih, k] h+t2) = J2^ Gxik, I] h) Gx{l, k; h). (2.32) 
I 

The measure factor I in the composition law comes from the circular nature of a 
strip. Writing the composition law for ti — 1 we see that the one step propagator 
acts as a transfer matrix, 

G,(/i, k-.t + l) = Y.l Gxik: I; 1) Ga(/, k; t). (2.33) 

In the following we derive Gx{k,k',t) by iterating (2.33) t times. The iteration 
procedure can be conveniently carried out by introducing the generating function 
for Gx{h,k;t) 

Gx{x,y;t) = J2^''y'Gxik,l;t), (2.34) 

k,l 

where x and y can be naturally interpreted as Boltzmann weights related to the 
boundary cosmological constants of individual triangles, 

X = e-^'", y = e-^°". (2.35) 

Analogously we write the Boltzmann weight related to the bulk cosmological con- 
stant as follows, 

g = e-^''\ (2.36) 

The above introduced notation implies that the total Boltzmann weight of one 
strip can be determined by associating a factor of gx with triangles that have the 
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spacelike edges on the entrance loop and a factor gy to triangles where the spacelike 
edges are on the exit loop. The one step propagator is now easily computed by 
standard generating function techniques as follows, 



G{x,y;g;l) = + ^ + ^ + ... 

k=l \m=l n=l / 

where the factor | comes from dividing by the volume of the automorphism group 
for periodic triangulations. Evaluating the summations in (2.37) we readily obtain 

n( ^\ f 9x gy Y , A gxgy \ , . 

From this expression it can be seen that the one step propagator with fixed bound- 
ary lengths is given by 

1 fh + h 



G{h,W,9;l) = 9'''-''j^^\^' 'J, (2.39) 

where the division by the volume of the automorphism group now makes its ap- 
pearance in the guise of the factor j^^- 

To compute the "finite time propagator" G{x,y; g;t), we rewrite the composition 
law (2.32) in terms of generating functions and obtain the following, 

/f dz dz' 1 
J2^i27i (1 _ ^(^' ^5 9; ti)G{z', y; g; ts). (2.40) 

By setting t2 = 1 and performing the contour integration over z we obtain 



G{x,y;g]t) = j 



dz' 



27ciz 



12 



-^Gxix,z; 1) 
dz 



Gx{z',y;t-1). (2.41) 

z=l/z' 



Inserting the expression for the one step propagator yields the desired iterative 
equation for G{x, y; g; t), 

G{x, y; g;t) = G (j^, y;g;t-l^- Gx {g, y;g;t-l). (2.42) 
The implicit solution of this equation can be written as 



32 



2 2D Causal Dynamical Triangulations 



where Ft is defined iteratively by 



The fixed point F as defined by Ft{x) = Ft_i{x) is given by 



, = (2.45) 

By well known methods one can use the fixed point to find the explicit solution to 
the iterative equation (2.44) 

where 

At = l-F2*+2, Bt^F-F^'+\ Ct^F^-F^'. (2.47) 

The complete finite time propagator is now obtained by substituting (2.46) in 
(2.43), yielding 

G{x, y; g;t)^- log (1 - Z{x, y; g; t)) , (2.48) 
where we have defined 

Z{x, y; g; t) = (l - ^) ^ ^^"'f ^ ^ (2.49) 

The region of convergence of this result as an expansion in powers of x, y, z is 

1^1 < I, < 1, \y\ < 1- (2.50) 
2.7.1 The continuum limit 

One of the central philosophies behind the method of dynamical triangulations is 
"universality". This concept is well known in statistical mechanics and plays a 
pivotal role in renormalization theory. As applied to the case at hand it means 
that the precise form of the discrete amphtude (2.48) should not be important for 
the physics of the system. The essential physics should for example not depend 
on the type of polytopes that are used to regularize the path integral. 
One of the prerequisites for such a scenario is that there exists a so-called "critical" 
hypersurface in the space of parameters of the regularized theory. Near such 
a region the theory exhibits correlations that are much larger than the size of 
the building blocks. If this happens the macroscopic physics is insensitive to the 
regularization and one can safely take the limit where the building blocks are 
infinitesimally small and obtain a more or less unique continuum theory. 
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Precisely how the above sketched scenario is rcaUzcd in two dimensional causal 
dynamical triangulations can be most easily understood by analyzing the one step 
propagator as a function of g, x, y. To identify the critical surface in the parameter 
space spanned by the three couplings g,x and y we compute the average number 
of triangles in one strip. This can be readily accomplished by taking the derivative 
with respect to g of the one step propagator since g can be seen as the fugacity 
for the number of triangles in the system, 

{Na) = 9^G{x, y;g;l) = l- + ^— -. (2.51) 

dg l-gx l-gy l-g{x + y) 

The critical region of the parameter space can now easily be identified by observ- 
ing that the total number of building blocks diverges if the couplings satisfy the 
following relation, 

\g{x + y)\ = l, (2.52) 

where it should be noted that besides this relation the couplings are also subject to 
the convergence condition (2.50). Near the critical region (2.52) the typical length 
scale of the system is indeed much larger than the length scale of the individual 
building blocks. To define the continuum limit of the theory we assume canonical 
scaling dimensions for the bulk and boundary cosmological constants, 

^ = ^,6-"'^, x = Xee-«^, y^ycB-'^''. (2.53) 

Since we have performed the Wick rotation one is dealing with real valued g, x, y. 
Therefore the natural critical values for these couplings are determined by (2.50) 
and (2.52) yielding 

gc^ 1/2, Xe = l, yc = l, (2.54) 
leading to the following scaling relations, 

g = ^e""'^, X = e-"^, y = e"'^^. (2.55) 

The continuum limit of the propagator can now be determined uniquely by insert- 
ing the scaling relations (2.55) into the regularized result (2.48), giving 

Ga{X, Y,T) = - log (1 - Z^X, Y, T)) , (2.56) 

where 

Za{X, F, T) = — — (2.57) 

(X + Va coth Vat) (y + Va coth Vat) 

Note that in the continuum limit we defined the normalization of the wavefunctions 
such that 

GAiX,Y,T) = aG{x,y;g;t). (2.58) 

The power of a is uniquely determined by the requirement that the continuum 
propagator satisfies a continuum analogue of the composition law (2.40). The 
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continuum expression for the propagator where the boundary lengths are fixed 
instead of the boundary cosmological constants can now be obtained by an inverse 
Laplace transformation of (2.56) with respect to both X and Y, 

where Ii{x) is a modified Bessel function of the first kind. One indication that 
this object is a good propagator is that it satisfies the expected composition law 

poo 

G'A(Li,L2;ri + r2)= / dLLGA{Li,L;T,)GA{L,L2,T2). (2.60) 

Jo 

From the propagator one can naturally define the time dependent disc function by 
shrinking one of the boundary loops to zero, 

Wa{L,T)^Ga{L,L'^0;T), (2.61) 

with 

Ga(L,L'^0;T)^ A^e-^^^°*^^^. (2.62) 

sinh VAT 

This object could be given a cosmological interpretation in that it is the probability 
amplitude to find a universe of size L that has existed for a time T. Upon doing 
an integration over time one obtains an object that has the natural interpretation 
of a Hartle Hawking wavefunction, 

dTWA{L,T)^^-. (2.63) 
Alternatively, this object is often called the disc function. 



2.7.2 Marking the causal propagator 

In most expositions on two dimensional causal dynamical triangulations the con- 
ventions differ from the ones used in the previous section. Often, one employs 
propagators that have a marked vertex on the initial boundary, 

G^^''\hM\^)^hGy{hM-l\ (2.64) 

where the superscript notation G^^'*^) denotes that the entrance loop of the propa- 
gator is marked, (7^°'^^ imphes that the exit loop is marked and G^^'^) means that 
both loops are marked. In principle all vertices of the triangulation are indistin- 
guishable, the introduction of a mark on the boundary implies that one vertex on 
the boundary is distinguishable from all others. The main virtue of the marking 
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is that it removes the measure factor in the composition law for the propagator 
(2.32) 

G'^^-''\h,k;h + t2) = °n^i,Mi) G^^''\l,k;h), (2.65) 

I 

which corresponds to the following composition law for the Laplace transformed 
propagator, 

G^;-''\x, y; ^1 + *2) = ^ 2^ "^x, z-'; t,)G^^-''\z, y; t,). (2.66) 

Many formulas simplify with this trick, but the price one pays is that the propaga- 
tor is no longer symmetric. Although the physical symmetry of the two boundary 
components of the causal propagator is not reflected in the marked expressions we 
discuss their properties for future convenience. 

Since we have already computed the one step propagator without a mark, we can 
find the marked version by simply taking a derivative with respect to the boundary 
cosmological constant, 

G^'-''\x,y,g,l)^x^G^;:'\x,y-g-l) = ^ (2.67) 

dx {1 - gx){l - g{x + y)) 

To find the marked propagator for finite t one can find an iterative equation anal- 
ogous to (2.42) , 

G^'-''\x,y;g;t) = ( -^,y; g;t - l) . (2.68) 

1-gx \i-gx ) 

Evidently, one can use this relation in a similar iteration procedure as described 
in section (2.7) to obtain G^^'^\x,y; g;t) and then apply the scahng relations to 
find the continuum propagator. Instead, one can also apply the scaling relations 
(2.55) directly to (2.68) and find a differential equation for the continuum marked 
propagator. 



dT 

where 



^Gr'{X,Y,T) = -^ W{X)g1-''\x,Y,T) , (2.69) 



W{X) = - A. (2.70) 

The notation W{X) is introduced to anticipate generalizations, since this differen- 
tial equation can also be solved for more general expressions than (2.70). Equation 
(2.69) has the form of a standard first order differential equation, so to solve it 
uniquely one needs to supply one boundary condition. The natural condition to 
impose on the propagator is that it reduces to the identity for the limit of zero 
proper time T, 

Gf°)(Li,L2;T = 0) = d{L, - L^), (2.71) 
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which has the following analogue in terms of boundary cosmological constants, 

G^l-''\X,Y;T = 0) = ^^. (2.72) 
Solving the differential equation (2.69) with this boundary condition gives 

^{i;0) (X YT)- ^^^^'^^^ ^ (n yoN 

where X{T; X) is the solution to the characteristic equation 

^^-W{X{T)), X{T^O)^X. (2.74) 

Since here we are interested in pure CDT, the definite form of this equation is 
known since W{X) — X^ — A and one can solve the characteristic equation explic- 
itly, 

XiT) = + (2.75) 

Inserting this expression in (2.73) it is easily verified that this form of the propa- 
gator is fully compatible with the previously derived results. If one marks the 
unmarked propagator (2.56) by taking a derivative with respect to the initial 
boundary cosmology constant X, the expression coincides precisely with the above 
derived result (2.73) . 

2.7.3 Hamiltonians in causal quantum gravity 

An interesting observation regarding the differential equation for the causal prop- 
agator (2.69) is that it can be viewed as a Wick rotated Schrodinger equation, 

-^G^;^'\X,Y,T) = HxG^^'\x,Y,T), (2.76) 

where the Hx is the quantum effective Hamiltonian and is given by 

^x = (^'-A)^ + 2X (2.77) 

By inverse Laplace transforming the Schrodinger equation (2.76) we can find the 
Hamiltonian in the "position representation" or, more appropriately the length 
representation, 

^^arked ^ -L—— +AL. (2.78) 



9L2 

This operator is selfadjoint on the Hilbert space T-l = C^{W^, L~^dL). Care should 
be taken however, since the above defined Hamiltonians are not derived from a 
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symmetric propagator. The proper Hamiltonian can be obtained by taking the 
continuum hmit of the nonmarked iteration equation (2.42), yielding 

which is selfadjoint on the Hilbert space H = £^(M+,L(iL), where the measure 
originates from the basic composition law (2.60). Physically, the Hamiltonian is 
well defined as it is bounded from below and it is even possible to find its spectrum 
explicitly. Solving the eigenvalue equation 

HlML) = Er^ML), (2.80) 

yields the eigenfunctions 

ML) = 2y^e-^^L^(2^/AL), (2.81) 

where Ll^{x) is a generalized Laguerre polynomial. Further, the eigenvalues are 
given by 

= 2VA(n+ 1). (2.82) 

Observe that the spectrum is equidistant, making the quantum mechanics of the 
problem rather similar to that of the harmonic oscillator. The main difference 
with the harmonic oscillator is that the Hilbert space of the CDT Hamiltonian is 
defined by the square integrable functions on the real half line L e [0, oo) and not 
on the whole real line L e (—00,00). 

To make contact with the continuum treatment of causal quantum gravity by 
Nakayama [30] we need to absorb the measure factor L in the normalization of the 
wavefunctions. The marking procedure described in section (2.7.2) does precisely 
this but the measure is absorbed in the initial state only, leading to a nonsymmetric 
propagator. In the continuum theory it is possible to absorb the measure factor 
in both the initial and final states by the rescaling ipn{L) = -^(pn{L). Hence, the 
Hamiltonian with flat measure is defined by 

H£^'iPn{L)=HL^ML), (2.83) 

giving 

m^"* = -L^ - -^ + — +AL. (2.84) 

To see how this Hamiltonian appears in the continuum derivation we consider the 
non-local "induced" action of 2d quantum gravity, first introduced by Polyakov [7] 



S[9] 



J dtdx^ Q-R^J-R^ + , (2.85) 
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where R is the scalar curvature corresponding to the metric g, t denotes "time" 
and X the "spatial" coordinate. 

Nakayama [30] analyzed the action (2.85) in proper time gauge assuming the man- 
ifold has the topology of a cylinder with a foliation in proper time t, i.e. the metric 
was assumed to be of the form: 

9 -(I J^X (2-86) 



^{t,x] 

It was shown that in this gauge the classical dynamics is described entirely by the 
following one-dimensional action: 

where 

l{t) = i y dx^, (2.88) 

and where k is an integration constant coming from solving for the energy-momentum 
tensor component Tqi = and inserting the solution in (2.85). 
Thus Leant = T^l{t) IS precisely the length of the spatial curve corresponding to 
a constant value of calculated in the metric (2.86). Nakayama quantized the 
actions for k — {m -\- 1)^, m a non-negative integer, and argued that in the 
quantum theory k = (m + |)^. The classical Hamiltonian corresponding to the 
proper time gauge action (2.87) was derived and reads as follows 

H^ = Iiimi + {m+\fj + M, (2.89) 



where is the canonical momentum conjugate. Subsequently, the quantum 

'dl 



Hamiltonian is found by the straightforward replacement — > — 



= -/^ - I + (m + i)' y + A L (2.90) 

Rather remarkably, upon the identification I ^ L we see that for m = the result 
coincides precisely with the flat measure Hamiltonian derived by CDT (2.84). An 
interpretation for the higher m quantum numbers in the context of CDT can be 
found in [55]. However, a complete understanding of the significance of these 
quantum numbers is in the opinion of the author still lacking. 
Recall that the parameter / in Nakayama's Hamiltonian is not the physical length, 
Lcont) it differs from L^ont by a factor tt. Consequently, it is natural that also in 
CDT the physical length is defined as -Lcont — '^L. This observation will turn out 
to be important for the emergent geometry discussion in chapter 4. 
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Despite recent progress [26, 27], little is known about the ultimate configuration 
space of quantum gravity on which its nonperturbative dynamics takes place. This 
makes it difficult to decide which (auxiliary) configuration space to choose as start- 
ing point for a quantization. In the context of a path integral quantization of grav- 
ity, the relevant question is which class of geometries one should be integrating over 
in the first place. Setting aside the formidable difficulties in "doing the integral" , 
there is a subtle balance between including too many geometries - such that the 
integral will simply fail to exist (nonperturbatively) in any meaningful way, even 
after renormalization - and inchiding too few geometries, with the danger of not 
capturing a physically relevant part of the configuration space. 
A time-honored part of this discussion is the question of whether a sum over 
different topologies should be included in the gravitational path integral. The 
absence to date of a viable theory of quantum gravity in four dimensions has not 
hindered speculation on the potential physical significance of processes involving 
topology change (for reviews, see [56, 57]). Because such processes necessarily 
violate causality, they are usually considered in a Euclidean setting where the 
issue does not arise. In our models for topology change we do wish to capture 
some Lorentzian aspects though. In section (3.2.1) we therefore address some of 
the causality issues and argue that the causal structure of the manifolds is only 
modified mildly. 

The focus of this chapter is devoted to the introduction of spatial topology change 
in two dimensional quantum gravity. This might be viewed as a bit of an ironical 
undertaking since the characteristic that makes the quantum geometry of CDT 
better behaved than its Euclidean rival is precisely its fixed spatial topology! We 
do however show that when a coupling constant is introduced the effect of the 
spatial topology changes is much less severe. A natural scaling for this coupling 
constant is presented, where the dynamics of an individual spacetime and splitting 
into baby universes both contribute to the continuum limit of the theory. Conse- 
quently, time scales canonically and the limit where the coupling constant tends to 
zero is smooth and continuous and gives back the results of the bare CDT model. 
This should be contrasted with Euclidean quantum gravity defined through dy- 
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namical triangulations, since in that case the continuum dynamics of the theory is 
completely dominated by the baby universes, making the geometry highly fractal. 
Rather than going into the details of the theory with the coupling constant straight- 
away we will first explain the role that baby universes play in the relation between 
Euclidean and Causal dynamical triangulations in section (3.1). The formalism 
of CDT can easily be extended to allow for spatial topology change whereby we 
admit the geometry to split into baby universes. It turns out that if one allows 
for the baby universes the splitting process will dominate the path integral com- 
pletely, wiping out all traces of the dynamics of each individual spatial universe 
[21]. At the same time it is shown that this continuum limit corresponds uniquely 
to the dynamics of Euclidean quantum gravity. Explicitly it is shown how to red- 
erive the Hartle Hawking wavefunction and propagator of Euclidean dynamical 
triangulations in the continuum limit. In this discussion we closely follow [58]. 
In sections (3. 2), (3. 3) and (3.4) we describe our published work where we generalize 
the above described construction [59]. Particularly, we associate a coupling con- 
stant that is reminiscent of the string coupling constant with the spatial topology 
fluctuations. We show that in a suitable scaling limit the spatial topology changes 
contribute to the path integral in a controlled manner without dominating the 
quantum geometry. 

3.1 Euclidean results with causal methods 

In section (3.1.1) we allow for spatial topology change and show how to explicitly 
introduce the baby universes in the discrete formalism and we compute the Hartle 
Hawking wavefunction in the continuum limit. The result is the well known disc 
amplitude of two dimensional Euclidean quantum gravity. 

3.1.1 Spatial topology change 

We now address the implementation of spatial topology changes by generalizing 
the discrete CDT framework. Although a multitude of constructions that realize 
topology change exist, they are all equivalent in the continuum limit as has been 
checked for some cases by the authors of [21]. In the following we introduce one 
specific realization and show that in the continuum limit the familiar results from 
Euclidean quantum gravity are recovered. 

The first step in the construction is to generalize the one step propagator, or 
transfer matrix, of CDT by alleviating the constraint on the spatial topology of 
its initial loop. Particularly, we include strips for which the entrance loop, say 
of length Zi, has the topology of a "figure eight". A natural procedure to create 
such a figure eight is to non-locally identify two points of a spatial universe with 
topology of an S^. Incorporating this pinching process leads to a factor of h in the 
combinatorics for the one step propagator since the pinching is allowed to happen 
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Figure 3.1: Illustration of a one step propagator with a "baby universe". 



at any of the h vertices. A baby universe is now created by associating one of 
the loops of the figure eight with the boundary of a disc function whilst the other 
loop is associated with the initial loop of a regular one step propagator (fig. 3.1). 
Combining all of the above, we can write the new transfer matrix in terms of the 
old, or "bare" , transfer matrix and the as yet undetermined disc function 



h-i 



Ga(/i, k; 1) = (/i, k; 1) + V kw{h-l, g) G'^'il, k; 1) 



1=1 



(3.1) 



where it should be noted that we have dropped the superscript notation for the 
marking and have used the notation Gx{li,k'A) = G^x {h, k'A)- 
This new, or "dressed" , one step propagator satisfies the same composition law as 
the bare transfer matrix (2.65), i.e. 



G\{li, k] ti + ^2) — Gx{li, I; ti)Gx{l, k', h)- 



(3.2) 



Particularly this implies that the dressed one step propagator can still be used as 
a transfer matrix despite the fact that the disc function contains an infinite sum 
over time steps 

Gxih, k; t) = J2 ^^(^1' 1) ^^(^' ^-1)- (3-3) 
I 

Performing a (discrete) Laplace transformation of eq. (3.3) leads to 



G{x,y;g;t) 
dz 



2m z 



Gf\x, z-'; l)+x-^ (w{x; g)Gf\x, z''; 1)) 



G{z,y;g;t-1), 



(3.4) 



or, using the explicit form of the transfer matrix G(^h){x, z; g; 1), formula (2.67), 

G{x,y;g;t)= 1 + x — \-xw{x,g)— G ,y;g]t-l). (3.5) 

L ox oxi 1 — ox VI — qx / 
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disc function 




disc function 



propagator 




Figure 3.2: Decomposition of the marked CDT disc function in another CDT disc 
function and a propagator. 



Note that this is not a closed equation, since so far neither the disc amphtude 
w{x,g) nor G{x,y] g]t) are known. Although this means that we cannot derive 
the discrete expressions, it will be shown by using scaling arguments that one 
can uniquely determine the continuum disc amplitude W\{X) and propagator 
Ga{X, Y;T). As in the case for CDT without topology change we assume canonical 
scaling for both the boundaries and the cosmo logical constant. 



9 



X 



-aX 



y 



-aY 



(3.6) 



In the following arguments no specific choice for the scaling of the time variable 
will be assumed as its scaling will be determined at a later stage. Interestingly, 
the composition law and the canonical scaling of the boundary lengths are enough 
to determine the scaling of the dressed propagator 



Gx{li, I2, t) — ^ a Ga{Li, L2; T). 



(3.7) 



Changing the boundary conditions from fixed boundary length to fixed boundary 
cosmological constant amounts to taking a Laplace transformation and implement- 



ing the canonical scaling of the boundary cosmological constant x 
to 



-aX 



leading 



G,{x,l2,t)'^GA{X,L,,T) 



(3.8) 



and 



Gx{x,y;t)^-^a-'GA{X,Y;T). 



(3.9) 



The scaling relation for the disc function is not as simple to obtain as the scaling 
for the propagator and we show that it actually depends on the scaling of time. 
Where for the propagator we use the composition law to derive its scaling relation 
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Figure 3.3: Decomposition of the disc function including spatial topology change. 



we use the following exact combinatorial identity to determine the scaling relation 
for the disc function 

9 ^^^^ = E E ^(^' ^' 9; t) I w{l, g), (3.10) 
9 t I 

or, after the usual Laplace transformation, 

dg ^ J 27riz ^ ' ' dz ^ ' 

These identities reflect the fact that if one introduces a mark anywhere in the 
bulk by taking a derivative with respect to g, the disc function can be decomposed 
into a propagator and another disc amplitude. The situation for the bare model 
is illustrated in fig. 3.2. In the figure on the right we have highlighted the points 
that have a distance t to the entrance loop. Since the bare model does not allow 
for spatial topology change these points constitute one spatial universe with the 
topology of an 5*^. In fig. 3.3 the situation is illustrated in the case where one 
does allow for topology change, in this case all points with equal distance t from 
the entrance loop form a spatial section with the topology of several S^'s and the 
mark is located on one of them. 

Since our current objective is to find the scaling for the disc function we assume 
the following general scaling ansatz 

w{x, g) = Wns{x, g) + aPWp^^X) + less singular terms. (3-12) 



In the case ry < the first term is irrelevant in the continuum limit, it does not 
appear in the computation of any continuum quantities. However, if r/ > a term 
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Wns will generically be present [60]. Additionally, the general ansatz for the scaling 
of the time variable reads as follows 

T = a't, £ > 0. (3.13) 

As shown in section (2.7.1) both space and time scale canonically in the bare 
model which corresponds to e = 1. Below we show that allowing the branching 
into baby universes to contribute to the continuum limit forces the scaling of time 
to be anomalous, creating an inherent asymmetry between the time- and space 
directions. 

Inserting the scahng relations (3.12) and (2.55) into eq. (3.11) we obtain 
dg OA ~ 



^J^tJ dZGj,{X,-Z;T) 



(3.14) 



where (x, (7) = {xciQc) in the non-singular part. 

Prom eq. (3.14) and the requirement e > it follows that the only consistent 
choices for -q are 

Scaling 1: ry < 

As can be seen from (3.12), this range of values corresponds to the situation 
where the non scaling part of the disc function is irrelevant and the physics 
is completely independent of the cutoff. 



dk 2^ 



= ^—f4Tf,lZ G,(X, -Z; T) (3.15) 

J J Zc oZ 



The continuum limit can be taken for any < since (3.15) does not depend 
on its explicit value. The value of e on the other hand is fixed and one needs 
to have e = 1 for the continuum limit to exist. Summarizing, if the scaling 
of the disc function is such that non scaling contributions are negligible 
in the continuum limit, the time variable automatically scales canonically. 
Obviously the bare CDT model falls in this class of scahngs since it has 
T] — —1 and e — 1. 

Scaling 2: 1 < rj < 2. 

For this class of scalings for the disc function formula (3.14) splits into two 
equations 



OA 2a= dz 



dT / dZ Ga{X,-Z;T), (3.16) 
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and 



dWns 



dg 



._f^!l/dT/.ZG.(X.-Z;T)i™ (3.17) 

9=9c J J c 

From (3.16) it follows that dP~'^ = and from (3.17) one sees that = 1. 
Combining these requirements we are led to the conclusion that e — 1/2 and 
T] = 3/2, which are precisely the values found in Euclidean 2d gravity. Let 
us further remark that eq. (3.16) in this case becomes 

dWA{X) ^ ^^^^^ ^^^^^ ^ ^g^^g^ 



dK 

If one rescales the couplings it is possible in general to absorb the constant 
originating from the non universal terms. This implies that the continuum 
dynamics of the theory strongly depends on the fact that there are non 
universal terms surviving the continuum limit but their precise value does 
not play a pertinent role. It should be noted that (3.16) expresses a different 
relation between the disc function and the propagator than was used for 
the bare model, it differs from (2.61) by a derivative with respect to the 
cosmological constant. Finally, inserting e = 1/2 and r] = 3/2 into eq. (3.17) 
yields 

J dT J dZ Ga{X, -Z-T) ^^q^^ = const, (3.19) 

where as in (3.18) the constant originates from the non scaling terms of the 
disc function and its value does not play a significant role in the continuum 
theory. 

The relation (3.18) possesses a remarkable interpretation in terms of baby uni- 
verses. Basically it states that near any mark in the bulk of the continuum Hartle 
Hawking wave function there is a baby universe at the scale of the cutoff. Since the 
location of the mark is arbitrary it implies that near every point there is a baby 
universe at the cutoff scale which is illustrated in fig. 3.4. This does not mean 
however that all baby universes are of negligible size. An additional indication 
that cutoff size geometries play an important role in the dynamics comes from the 
Laplace transform of (3.19) 

J dT J dL Ga{L,L';T) ll Wa(L') = const, x 5(L), (3.20) 

which shows that the distribution of geometries is such that it is strongly peaked 
around universes that have minimal boundary length. In the following we show 
that these microscopic artifacts are an important feature of 2d Euclidean quantum 
gravity since the CDT model with baby universes exactly reproduces the contin- 
uum equations of the Euclidean model in the scaling limit. Furthermore, to make 
contact with the Euclidean theory, (3.18) is used in an essential way. 
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Figure 3.4: At every point in the quantum geometry there is an infinitesimal baby 
universe. 



Having derived the scahng relations we can now analyze the scaling limit of (3.5) 
and find an equation for the dynamics of the propagator. In order for the equation 
to have a scaling limit at all, Xc, Qc and Wns{xc,gc) must satisfy two relations 
which can be determined straightforwardly from (3.5). The remaining continuum 
equation reads 



df 



Ga{X,Y;T) 



—a 



_d_ 
dX 



d 



dX 



(X^ - A)GAiX,Y;T) 
Wa{X)Ga{X,Y;T) 



(3.21) 



The first term on the right-hand side of eq. (3.21) is precisely the one we have 
already encountered in the bare model, while the second term is due to the creation 
of baby universes. In case where r] < 1 the first term in (3.21) is subdominant and 
the scaling will not be compatible with baby universes in the continuum limit. So 
from (3.21) it can be seen that for 77 < 1 we need to have that e = 1 leaving us 
with the continuum differential equation for the propagator of the bare model 



^GAiX,Y;T)^-^ 



(X2-A)Ga(X,F;T) 



(3.22) 



where we remind the reader that this differential equation can naturally be inter- 
preted as a Wick rotated Schrodinger equation of a single string propagating with 
respect to its own time on world sheet. 



-^GAiX,Y;T) = HxGAiX,Y;T). 



(3.23) 
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Prom the Laplace transform of this equation one sees that if is a simple Hamilto- 
nian containing a kinetic term, a potential induced by the cosmological constant 
and no interaction terms, 

Jjn^arked ^ -L—^ + A L. (3.24) 

For the second scaling however, where 1 < < 2, the last term on the right-hand 
side of (3.21) will always dominate over the first term. Consequently, the first 
term does not survive the continuum limit leaving one with an equation without 
a Hamiltonian containing a kinetic term. So for this scaling the dynamics of the 
world sheet is governed purely by the interactions of splitting strings. Equivalently, 
we can say that once we allow for the creation of baby universes, this process will 
completely dominate the continuum limit. As we have seen from (3.16) and (3.17), 
(?7, e) = (3/2,1/2) is the only consistent scaling that allows for baby universes. 
Inserting this scaling in (3.21) we obtain the following continuum equation 

^ Ga(X, Y; T) = - A ^Wa{X)Gj,{X, Y; T)] , (3.25) 

which, combined with eq. (3.18), determines the continuum disc amplitude W\{X). 
Integrating (3.25) with respect to T and using that Ga(-^i, -^2; 2^=0) = 5(Li— L2), 
i.e. 

G'a(X,L2 = 0;T = 0) = 1, (3.26) 

we obtain 

Prom dimensional analysis one can easily see that Wl{X) — X^F{y/A./X). This 
imphes that the solution of the disc function reads as follows 

Wa{X) = V-2AA: + 62^3 + c2A3/2. (3.28) 

However, not all values for b and c are physically acceptable. The inverse Laplace 
transform of (3.28), Wa{L), should be bounded for all L > 0. This constraint gives 
the following expression for the disc function 



where the constant 6 is a constant that reflects specific details of the discrete 
statistical model, as described by equation (3.18). We discussed above that this 
constant does not have an obvious physical significance since it can be absorbed 
into the cosmological constant. Absorbing the irrelevant constant b one obtains 
the disc function wj;^^\x) of 2d Euclidean quantum gravity. 



Wa{X) = {X - l-^/A)\/ X + Va. (3.30) 
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Figure 3.5: A representation of the global aspects of a trouser geometry with 
modular parameters {ti,T2,T3). 



Note that the defining equation for the disc function (3.25), can alternatively be 
derived by several methods within 2d Euclidean quantum gravity [37, 60, 61]. In 
those computations it is clear that T can be interpreted as the geodesic distance 
between the initial and final loop. 

For a more detailed account on baby universes in two dimensional Euclidean quan- 
tum gravity the reader is referred to [62] . 

3.2 Introducing the new coupling constant 

In this section we present a model which can be viewed as a hybrid between 
Euclidean and causal quantum gravity. Recall that the pure CDT model possesses 
a Hamiltonian that governs the dynamics of an individual string, but that there is 
no interaction term in the dynamical equation to allow for spatial topology change 
(3.22). In the case of the Euchdean theory the opposite is true, the equation 
governing the dynamics of the propagator only contains an interaction term and 
no kinetic or potential terms for an individual string (3.25). In the previous section 
it was shown that both theories can be obtained by two different scaling limits of 
one unifying statistical model based on CDT (3.1). The current objective is to 
show that there is a natural adaptation of this unifying model such that its scaling 
limit acquires single string dynamics while at the same time allowing interactions 
in the form of spatial topology change. The essential ingredient in accomplishing 
this is to introduce a new coupling in the statistical model (3.1) whose scaling is 
such that both the interaction term and the dynamical term in (3.1) contribute in 
the continuum limit. Before discussing the construction of the model we examine 
some aspects of geometries with spatial topology change. The simplest example of 
a two dimensional geometry with spatial topology change is the so-called "trouser" 
geometry. A trouser geometry is a manifold with three boundary components with 
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figure eight 




T 



Figure 3.6: A representation of the global geometry of a trouser geometry in causal 
dynamical triangulations 



the topology of an S^. In general the so-called legs of the trousers can have any 
length which can be parameterized by three modular parameters as illustrated in 



We are interested however in a more restricted subclass of trouser geometries which 
are of the form depicted in fig. 3.6. We demand that we have an initial boundary, 
where every point on this boundary has the same distance to the two final loops. 
This restriction means that on top of the length of the boundary components we 
need just two parameters to describe the global characteristics of the geometries, 
one describing the length of the trouser leg belonging to the initial boundary 
component and one that specifies the length of the trouser legs related to the final 
boundary components, since we required the final boundary components to have 
the same distance to the initial boundary. 

Now we ask ourselves the question, what is the contribution of such a geometry 
to the path integral? As described in chapter 2 the action for two dimensional 
gravity contains just a volume term proportional to the cosmological constant and 
a topological term coming from the two dimensional Einstein-Hilbert action and 
the Gibbons-Hawking- York boundary term. 



where xi^) = (2 — 2g — 6) is the Euler characteristic of the manifold. The weight 
of an individual geometry in the path integral is therefore proportional to the expo- 
nential of the volume and the exponential of the number of handles and boundary 
components. 



We see that geometries with the topology of a cylinder are not weighted by gs, since 
the Euler characteristic of a cylinder is zero. The Euler characteristic of a trouser 



fig. 3.5 . 



Sm = AVm + ^X{M) 



(3.31) 




(3.32) 
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geometry is one however, since the trouser geometry has one more boundary com- 
ponent than the cyhnder. Consequently the weight of the trouser geometry in the 
path integral is not only determined by the cosmological constant, but also by 
the coupling gs- For the purposes of this thesis we are not particularly interested 
in computing amplitudes with the topology of a trouser, but we are interested 
in calculating propagators, i.e. amplitudes that have two macroscopic boundary 
components. If one shrinks one of the final boundary components of a trouser 
geometry to zero one obtains a geometry with two boundary components that we 
interpret as a "cylinder geometry with one baby universe". Shrinking a bound- 
ary component does not alter the Euler characteristic, allowing us to conclude 
that a geometry with n baby universes contributes to the path integral with a gg 
topological weight. Note that in the derivation of the disc amplitude of Euclidean 
dynamical triangulations of section (3.1) we did not associate a coupling constant 
to the baby universes. Subsequently, the formation of baby universe was not sup- 
pressed with the result that the continuum limit is dominated by cutoff scale baby 
universes. In section (3.3) we show that adding the coupling constant makes the 
term baby universes a bit misleading since the outgrowths have a certain finite size 
in the continuum limit and can be made arbitrarily big or small depending on the 
value of gs- 

3.2.1 Lorentzian aspects 

The reasoning in the previous section was based on a picture where the metrics 
are Wick rotated from Lorentzian to Euclidean signature. Performing an inverse 
Wick rotation is not so easy however for geometries with baby universes, since 
these geometries do not admit Lorentzian metric everywhere. The reason for this 
comes from the fact that it is not possible to find a non vanishing vector field 
everywhere. It is however possible to find a vector field everywhere except for 
a finite number of points. For some values of the time parameter the universe 
develops a baby universe and the topology of a spatial universe splits up into two 
different components. Precisely at the moment of a split the spatial geometry has 
the topology of a figure eight and we see that the central point of the figure eight 
is a saddle point if one views the geometry as being embedded in M"^. This point 
is called a Morse point in the mathematically oriented literature [57]. Contrary 
to a generic point on the manifold, such a Morse point does not possess a unique 
timelike vector field perpendicular to its spatial slice. In fact it has two future 
directed normal timelike vectors, showing that the neighborhood of such a point 
has an anomalous causal structure. More precisely, such a point has two future, 
and two past light cones which we refer to as a double light cone structure [57]. 
Consequently, if one universe splits in two, each of the resulting universes carries 
a light cone belonging to the Morse point as is illustrated in fig. 3.7. Conversely, 
the universe carries the two past light cones of the Morse point before the split. 
Because of the peculiar causal structure around the Morse points we are not quite 
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Figure 3.7: Illustration of the double light cone causal structure around a Morse 
point. 



sure whether the usual definition of the Wick rotation in CDT is also legitimate 
here. There are some indications that the Einstein Hilbert action develops complex 
valued singularities at these points [63]. Although the Wick rotation at these points 
presents one with an interesting conundrum, we confine ourselves to the Euclidean 
sector for the remainder of this thesis. 



3.3 Dynamics to all orders in the coupling 

In this section we discuss the explicit construction of the model with a coupling 
constant for the spatial topology changes and show that it can be solved to all 
orders in the coupling constant in a suitable continuum limit. 
Our starting point is equation (3.1) with the addition of a coupling to the inter- 
action term, where in section (3.2) we argued that the coupling can naturally be 
denoted as gs, 

h-i 

Gxih, k; 1) = (Zi, k; 1) + 2g, ^ l,wx,gAh-l) G'J)(/, k; 1). (3.33) 

1=1 

After a discrete Laplace transform one obtains 

G{x,y;g;t) = 



G{z,y]g]t-\). 

(3.34) 



If the baby universes are to survive the continuum limit the coupling constant in 
a controlled manner we have to scale gs with the lattice cutoff a. To find the 
appropriate scaling for the coupling we implement the following ansatz 

gs = a^gs + less singular terms. (3.35) 
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Using the explicit form of the transfer matrix G 



X, z 



-1. 



1) and inserting the 



scahng relations (2.55), (3.35) into expression (3.34) one is led to the following 
continuum equation, 



A)Ga{X,Y;T) 



dX 



WA,g,{X)GA{X,Y-T) . (3.36) 



If we demand that the first term on the right hand side of eq. (3.36) survives the 
continuum limit we need e = 1 as in the bare CDT model. Contrary to (3.21) 
we have the additional freedom to scale the coupling constant gs- This enables 
us to adjust the scaling so that also the second term survives the continuum limit 
yielding rj + C, = 2. Inserting this relation in (3.36) wc obtain the desired dynamical 
equation that contains both the dynamical term for propagation of single strings 
and a term governing string interactions. 



^G^'''\X,Y,T) 



d 



dX 



[{{X' - A) + 2gsWA,,,{X)) G^'-''\X, Y, T)] . (3.37) 



We observe that this equation is of the same form as the equations for the CDT 
models without a coupling constant (2.69), 



d 

df 



0^m(^X,Y;T) 



d 



dX 



Wa,,s{X)G^'-''\X,Y;T) 



but with a different form of W{X), 

Wa,,,{X) = {X' - A) + 2gsWA,,,{X). 



(3.38) 



(3.39) 



Note that requiring both terms in (3.37) to survive the scaling limit does not fix 
the scaling uniquely. However, since we merely introduce a coupling constant to 
the baby universe model and did not introduce any new configurations in the path 
integral, the disc function should still satisfy (3.10). As before we can use this 
relation to further constrain the scaling. The dynamics of our new model requires 
£ = 1, so we are led to the conclusion that the model falls into the first of the two 
scaling classes defined in section (3.1) implying 77 < 0. Consequently, the relation 
between the continuum disc function and the continuum propagator is the same 
as for the bare CDT model. 



dWA,,,{X) 
OA 



^ J dT J dZG^^'^\X,-Z;T) 



dW, 



A,9s 



dZ 



(3.40) 



It seems that we are still left with a predicament since we are unable to specify 
the exact scaling for the disc function. In fact it is only an illusory problem, since 
the disc function only appears in the dynamical equations in combination with 
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the coupling constant as gsW. Therefore, we in principle only need to determine 
the scaling of this combination which means that the relevant part of the scaling 
is already determined by the dynamical equation itself. In the model we consider 
however, the scaling is completely fixed by requiring the disc function to reduce 
to the disc function of the bare CDT model for zero gs, implying 

'UJx,gsix) =(i~^^A,9si^)^ 9s = a^gs- (3.41) 

In the remainder of this section (3.40) and (3.37) are used to derive a differential 
equation for WA,gs{X). The equation will be rather implicit however, since the 
equation does not only depend on W\^gg{X), but it also depends explicitly on 
the solution of the equation W\^gg{X) = 0. Remarkably, one is able to solve the 
equation uniquely provided the disc functions Wa.^s {L) satisfy the natural physical 
requirement that they fall of at infinity. As a first step we solve the dynamical 
equation (3.37) to obtain the propagator in terms of W\^gg{X), 

G('-''\X,Y,T) = ^^:^siXiT)) 1 _ 

Wj,,g,{X) X{T)+Y ^ > 

Inserting this into the consistency condition (3.40) and performing the integration 
over Z we obtain 

dW^,g,{X) _ r WA,,siXiT)) dWA,g, 

To evaluate the integral we conveniently use the characteristic equation and convert 
the integral over time into an integral over X, 

^^-Wj,,g,{X{T)). (3.44) 
Applying the characteristic equation to (3.43) we obtain 

= -^J— 7«^, (3.45) 

where one can easily evaluate the integral since the integrand is a total derivative, 

dWA,g,(X) _ WA,g,{X) - WA,g,(X^) 



(3.46) 



This equation can be rewritten as an equation for W\^gg{X) by reexpressing the 
disc functions as 
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giving 



or equivalently, 



(3.48) 



dA 



2{X^-X'). (3.49) 



This is the defining equation for W\^g^,{X) that we intended to derive. The simple 
form of the equation is deceptive as X^ is in fact a completely unknown function of 
gs and A, the only thing we know is that it is defined as the solution of WA^gg{X) ~ 
0. Another unknown arises when we try to solve (3.49) by integrating both sides 
with respect to A, the integration "constant", C{gs,X), can in principle be a 
general function of gs and X. So we conclude that (3.49) only determines W\^gg{X) 
up to two functions. Explicitly, the integral of (3.49) can be written as follows, 

W^,gAXr = AV {^) - 2X'A + X'h (^) , (3.50) 

where we have used dimensional analysis to parameterize the integration constant 
and the integral of X^ by two dimensionless functions, 

AV (-|,) = 2 p dA!X^ (gs, ^/A^') , C{gs, X) = X% (^) . (3.51) 
Below we use the following series expansions 

oo oo 

/te)=E^(;&)^ hm=Y.Km\ (3.52) 

n=0 n=0 

We show that the lowest order coeflicients of the expansions are determined by 
requiring consistency with the bare CDT model. Amazingly, all higher order coef- 
ficients are uniquely determined by demanding that the inverse Laplace transform 
of the disc function W{X) falls of at infinity at arbitrary order in the gs expan- 
sion. From the definition of W{X) in the pure model (2.70) one readily sees that 
to lowest order in gs we need to have 

WK,g,{Xf ^X^- 2X^A + A^ + 0{gs), (3.53) 

which implies /iq = 1, /o = 1 when comparing to (3.50). Using this result we write 
the disc function as follows. 



-{X' - A) + {X' - A)p + ^ ^ 

WA,gs{X) = ^2^^ . (3.54) 
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If we expand this expression to lowest order in gs we obtain 

'yA..(X) = f'^,+/'/^ + Ote). (3.55) 

By equating this result to the marked disc function of the bare model 

W^,o{X) = — (3.56) 
X + V A 

we conclude that hi — A, fi — —4. To obtain the higher order coefficients we 
proceed by analyzing the gs expansion of the disc function order by order. The 
first order correction to the disc function is obtained by inserting hi — A, fi — —4 
into (3.54) and expanding to first order in the coupling, 

W.,s(X) = ^ + 9. (-^JJZS^fc^ + '■""'xtT'^'" ) + (3.57) 
After an inverse Laplace transform we obtain 

VaVa,i(L) = ie^^(/2+/i2-l)-^+^e-^^ (2AL' + 2^/AL -f,-h, + l), 

(3.58) 

where we have introduced the following notation 

oo 

W^,,,{X) = J2WA,nWs- (3.59) 

n=0 

Demanding the disc function to fall of at infinity implies that the terms propor- 
tional to L and must vanish, leading to /i2 = 0, /2 = 1. So the result for the 
disc function at first order in the coupling reads as follows. 



WaAL) = (3.60) 

which can be confirmed by explicitly computing the "Feynman diagram" where the 
spatial universe is allowed to spht once. Obtaining the higher order coefficients is 
a bit messy, but the iterative procedure to compute them is completely analogous 
to the calculation for /i2 and /2. At each order of the gs expansion the disc 
function has the same form as (3.60). Every W\^n{L) contains three terms, a 
polynomial term in L, a term proportional to and a term proportional to 

g-VAL Demanding the disc function to be bounded at infinity implies that the 
polynomial and 0{e^^) terms should vanish. If one additionally uses the known 
results for hn-i and one obtains the /i„ and /„ coefficients uniquely. It turns 
out that boundedness of the disc function is such a stringent constraint that we 
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need hn ~ for all n >2. Inserting the nonzero coefficients of h (^) into (3.50) 
we obtain 

WA,gsiXY = X^- 2X''A + AgsX + F{gs, A), (3.61) 

where 

Figs,A) = A'fiq), q=^. (3.62) 

As noted above, demanding the disc function to be bounded at infinity also fixes 
all /„. Contrary to the /i„ coefficients however, the /„ coefficients are rather 
nontrivial but can be computed by the algorithm sketched above. Employing a 
symbolic computer program such as Mathematica one can readily compute the 
first coefficients (~ 20). Making use of Sloan's database of integer sequences it is 
possible to find a closed analytical expression for the coefficients /„ in terms of 
Euler gamma functions, 

r(|n-2) 

/o = 1, = ^ n>l. 3.63) 

' r(in + i)r(n)' - ^ ^ 

Given these coefficients we can sum the Taylor expansion and obtain the full non 
perturbative result for / {q) 

m = l + hF^{-i -h h 573^) - 42i^ (~i h h 573^) • (3-64) 

This means that we have solved our model and we can present the disc function 
with a non perturbative sum over spatial topologies, 

-iX^-A) + ^X^- 2X-K + AqVtx + A^/ jq) 
WA,gs{X) = — . (3.65) 

3.4 Relation to random trees 

In this section we give an alternative derivation of the disc function dressed with 
spatial topology changes (3.65) to show robustness of the result and to give more 
insight into the details of the quantum geometry. Particularly, the derivation we 
present below highlights the random tree structure of the configurations. To make 
the connection as clear as possible we start by deriving the one point function 
for a random tree model. Sometimes this one point function for random trees is 
referred to as a partition function for rooted branched polymers. The random tree 
model we consider is a statistical model consisting of edges that arc weighted with 
a fugacity z and three-valent vertices with a coupling constant A. A convenient 
way to view the partition function w{z,g) is that it is the generating function for 
the number of random trees, 

oo 

w{z,X) = Z^W2m{>^z'^)"' , (3.66) 
m=0 
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z z 

Figure 3.8: Pictorial representation of the iterative equation for the one-point 
function w{X, z) for branched polymers. 



where we use the fact that every vertex except the initial, or marked, vertex is 
accompanied by two edges. An easy way to evaluate the partition function w{z, A) 
is to notice that it should solve the following equation 

w{z,\) ^ z + \zw{z,\f. (3.67) 

The interpretation of this equation is illustrated in fig. 3.8. Equation (3.67) has 
two solutions but only one solution is compatible with the initial condition that 
there is only one tree with one edge, i.e. W2m, = 1- This solution is given by 

A) = — , (3.68) 

which has the following series expansion, 

w{z, \) = zT -r^^ (A^^)"^ . (3.69) 
^ ^ ^ m!(m + 1)! ^ ' ^ ^ 

m=0 ^ ' 

So we see that the number of random trees with m edges, W2m; is given by the 
m}^ Catalan number. This random tree model can be regarded as a simple model 
possessing some features of 0^ theory. Indeed it is possible to formulate scalar 0^ 
theory in as a sum over connected diagrams in the same way as the random tree 
model. One merely needs to replace the fugacity of the edges z by the standard 
scalar propagators. Here we view the random trees not as a model for particles 
interacting in an ambient space, but as a simple model for a theory of "interacting 
universes" in two dimensional quantum gravity. The idea is similar to the basic 
idea behind string theory, we "blow up" the Feynman diagrams by replacing the 
propagators of an interacting field theory by the propagators of a string theory as 
is illustrated in fig. 3.9. In our case this means that we replace the rather trivial 
propagators of the branched polymer model z by propagators that we computed 
from CDT (2.73). We remind the reader that the resulting model can either be 
viewed as a toy model for quantum gravity with topology change or as a string 
theory without a target space. To make the analogy as close as possible we write 
the equation for the generating function of the branched polymer as follows. 
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Figure 3.9: Illustration of spatial topology change as a string like generalization of 
a particle interaction. 



w{z, A) = wo{z) + XGo{z)w{z, Xf, (3.70) 

where Wo{z) and Go(^) denote the one point function and the two point function 
respectively, to zeroth order in the coupling A. From (3.67) we observe that the 
expressions for these objects coincide wq{z) = Go{z) = z. This is an understand- 
able coincidence, since the branched polymer is a very simple model. The equation 
analogous to (3.70) for our model of interacting spatial universes based on CDT is 

W^,g,{L) = WA,oiL) +gsj dTdL,dL2Gl;i\L, L, + T)W^a,,,(Li)W^a,,,(L2). 

(3.71) 

This is the defining equation for the disc function W{L). It is similar to the way 
equation (3.70) defines the one point function of the branched polymers w{z,X), 
since the disc function is the CDT analogue of the branched polymer one point 
function w{z,X). Instead of solving (3.71) we work with its Laplace transformed 
analogue, 

WA,gAX) = WA,o{X)+gs J dZ,dZ2G^li{X;-Z,,-Z2)Wj,,g^{Z,)WA,gsiZ2), 

(3.72) 

where the propagator is given by 

r^{i,i)fX- Y Y-T) = ^^'^^^^ \ _ ^ _ + ^ 

w^A,o(x) Ux + ri)2(x + r2) (x + ri)(x + F2)2. ■ 

(3.73) 

Performing the integrations over Zi and Z2 in (3.72) gives 

W,,AX) - «V„(X) (3.74) 
As in (3.43) it is convenient to use the characteristic equation to convert the 



3.4 Relation to random trees 59 



integral over T into an integral over X{X, T), 



W,„AX) = W,AX) + -3-^ / (3.76) 



X 



Seeing that the integrand is a total derivative one easily obtains 

Wa,o{X) 

Since this is a second order polynomial equation for W{X) it is readily solved and 
we obtain, 

^gs 

where 

WA,gsiX) = y#A,o(X)2 + 4^5 (l^A,o(X)VrA,o(X) +^siyA,,s(^oo)2). (3.78) 

Note that this equation was derived independently of the particular form of W\^q{X) 
and Vf^A,o(-'^)- Currently, however we are interested in an interacting model based 
on CDT. So we require the disc function and the propagator to reduce to the 
results obtained in the bare CDT theory, 

Wa,o{X) ^ #A,oW=^'-A, X^ = ^. (3.79) 
X + VA 

Inserting into (3.78) gives, 

Wj,,,,{Xr = {X' - + 4gs ((X -VA)+ gsWA,,,{VAr ) , (3.80) 
which can be conveniently written as 

WA,gs{Xf ^X^- 2X^A + 4gsX + F{gs, A), (3.81) 

where 

F{gs, A) = A^ - igsVX + iglWA,gs{VAf ■ (3.82) 
Given (3.81) and (3.39) we can obtain the following result for the disc function, 

W,.,AX) = -'^^ - ^' + ^""^ - l^l^ ^ + (3.83) 

which is exactly of the same form as derived in the previous section (3.3) but 
so far we have not yet determined the precise form of F(gf5,A). To abbreviate 
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the notation, we can convert the equation to dimensionless units by dividing all 
dimensional quantities by appropriate powers of the cosmological constant 

^ = A' ^=*' ^ = ^^ (3-84) 

and obtain 

= -<-'-^)+'^»W , (3.85) 

with 

LOqix) = Vx4-2x2 + 4gx + /(g). (3.86) 

In the derivation of section (3.3) the disc function is expanded in powers of q and 
the Taylor coefficients of F{q) are uniquely determined by demanding that the 
disc functions decay at infinity at every order of the expansion. Below we take a 
different route to obtain the /(g) by again using the characteristic equation in an 
essential way. Surprisingly, we find f{q) in a form that appears very different from 
the results of section (3.3), but is in fact exactly the same. Let us recall that the 
characteristic equation (3.44) can be used to define the time variable in terms of 
u;q{x) by 

X 

dx 



(jj{x) 



(3.87) 



Notice that we can only integrate t all the way to infinity if ujq{x) has a simple zero 
since x^o is defined as the solution of iliq{x) = 0. Together with (3.86) it implies 
that u;q{x) — should be of the following form. 



Uq[X) = [X 



c)y/{x + c+){x + c-), (3.88) 



where c, c+ and c_ are all functions of q that we determine below by taking the 
square of (3.88) and equating it with the square of (3.86). This gives us four 
equations for the four unknown functions c(g), c_|_(g), c_(g) and f{q), one equation 
for each power of x, enabling one to solve the system completely. If one solves the 
equation belonging to x^ one can eliminate one function and we can write 

Uq{x) — {x — cy{x + c + \/u){x + c — V^)) (3.89) 

where c+ = c + ^/u and C- — c — y/u. If we now expand equation (3.89) in powers 
of X and than equate it to the square of (3.86) we obtain 

- {2c^ + u)x^ + 2cux + c^{c^ -u) =x^ - 2x^ + Aqx + f. (3.90) 

Prom this we extract three equations, two simple relations expressing / and u in 
terms of c, 

/ = c2(3c2 - 2), (3.91) 
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u^2- 2c^ (3.92) 
and a third order polynomial equation for c 

c^-c = -q. (3.93) 

Inserting the solution of the polynomial equation for c in (3.91) gives f{q) which 
together with (3.86) allows us to find the complete solution for i2i{x). The solution 
of (3.93) is most conveniently expressed in terms of g = q as follows, 



VSc{q) = Zq + -, (3.94) 

where 



Zq 



-q+^q^-l)\ (3.95) 
Combining (3.94), (3.95) and (3.91) we obtain. 



/(?) = ^ ha + + 2 + ^ + ^ ) , (3.96) 



which appears to be very different from the expression found in the previous sec- 
tion, 

m -l + hFi (-1, -§; i; - 42i^ (-i h 1; , (3-97) 

but when one compares the Taylor expansion of both expressions we see that they 
are fully equivalent. 



3.5 Summary 

To set the stage for our model we reviewed the known relation between Euclidean 
and causal quantum gravity defined by dynamical triangulations in section (3.1). 
Some results of Euclidean quantum gravity can be derived by generalizing the 
formalism of causal dynamical triangulations to allow for spatial topology change. 
No "energy penalty" is associated with these topological fluctuations, manifesting 
itself by the fact that infinitesimal baby universes dominate the path integral. The 
fractal nature of the quantum geometry is rcficctcd by the non canonical values 
for both the "time" variable and the Hausdorff dimension. 

In section (3.2) we introduced a coupling constant for the spatial topology changes. 
Prom the Einstein Hilbert action of an elementary manifold with a change of spa- 
tial topology, the trouser geometry, we argued the naturalness of such a coupling. 
Additionally, we examined the geometry around the point where the spatial topol- 
ogy change occurs, the Morse point, and recalled that the causal structure around 
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such a point is non standard. Particularly, the Morse point features a doubling of 
the light cone structure, it possesses two past and two future light cones. 
In section (3.3) we discussed the details of the construction our new model of two 
dimensional quantum gravity with spatial topology change. The discrete kine- 
matical structure of the model is similar to the introductory section (3.1), the 
continuum behavior on the other hand is completely different. The presence of 
the coupling allows us to define a continuum limit where manifolds with all spatial 
topologies contribute but complicated topologies are suppressed by powers of the 
coupling constant. Especially, we were able to derive the disc function to all orders 
in the coupling constant and sum the series uniquely! 

An alternative derivation of the disc function dressed with topology fluctuations 
was presented in section (3.4). We showed that the disc function of the model can 
be derived from an iterative equation that is very similar to the generating function 
equation that deflnes the one point function of rooted random trees, or equivalently 
branched polymers. Besides providing additional insight into the structure of the 
quantum geometry of the model we also found that the hypergeometric functions 
that appeared in the result of (3.3) can be written in terms of a solution of a third 
order polynomial equation. 



4 

Hyperbolic space 



In this chapter, which is based on [64], we go back to an analysis of pure CDT. 
As in most approaches to quantum gravity, the CDT method was originally de- 
veloped for the quantization of compact manifolds. In the following however we 
generalize the boundary conditions of two dimensional CDT so that the typical 
geometries in the path integral have an infinite volume. Particularly, it is shown 
that given such boundary conditions a classical geometry with constant negative 
curvature and superimposed quantum fluctuations emerges from the background 
independent path integral. Furthermore, one can choose the boundary conditions 
such that the relative fluctuations become small in a concrete manner. To the 
knowledge of the author this is one of the few cases where a semiclassical geom- 
etry emerges from a genuinely background setup that can at the same time be 
studied by analytical methods. Another example were a semiclassical background 
emerges dynamically from a background independent scheme is four dimensional 
quantum gravity defined by causal dynamical triangulations. This model is too 
complicated to study with analytical methods, nevertheless it is one of the most 
promising attempts to formulate a realistic theory of quantum gravity (for a recent 
account see [29]). 



4.1 Non compact manifolds 

As mentioned before, 2d quantum gravity is intimately related with the study 
of non-critical string theories. The studies where the quantum gravity aspect has 
been emphasized mostly consider two dimensional Euclidean quantum gravity with 
compact spacetime. The study of 2d Euclidean quantum gravity with non-compact 
spacetimc was initiated by the Zamolodchikovs (ZZ) [20] when they showed how 
to use conformal bootstrap and the cluster-decomposition properties to quantize 
Liouville theory on the pseudo-sphere (the Poincare disc). 

Martinec [10] and Seiberg et al. [65] showed how the work of ZZ fitted into the 
framework of non-critical string theory, where the ZZ-theory could be reinterpreted 
as special branes, now called ZZ-branes. Let PF^(X) be the ordinary disc ampli- 
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tude for 2d Euclidean gravity on a compact spacetime. X denotes the boundary 
cosmological constant of the disc and A the cosmological constant. It was found 
that the ZZ-brane of 2d Euchdean gravity was associated with the zero of 



W^{X) = (X - X + VI. (4.1) 

At first sight this is somewhat surprising since from a world-sheet point of view 
the disc is compact while the Poincare disc is non-compact. In [66] and [67] it was 
shown how it could be understood in terms of world sheet geometry, i.e. from a 
2d quantum gravity point of view. When the boundary cosmological constant X 
reaches the value X — a/A/2 where the disc amphtude Wj;^{X) — 0, the geodesic 
distance from a generic point on the disc to the boundary diverges, in this way 
effectively creating a non-compact spacetime. 

Here we show that the same phenomenon occurs in two dimensional quantum 
gravity from causal dynamical triangulations. 



4.2 The hyperbolic plane from CDT 

Recall that the propagator with one mark on the initial boundary is given by 

G.(x,r;r)^aZ^^^^lj^, (4.2) 

where X{T, X) is the solution of the characteristic equation 

^ = -(X2-A), X{0,X)^X, (4.3) 

giving 

X{t,X) ^VXcoth^/A{t + to), X = VXcoth v^^o- (4.4) 

Note that although different in appearance this expression is equivalent to (2.75). 
Viewing Ga{X,Y;T) as a propagator with X and Y as coupling constants, X{T) 
can be viewed as a "running" boundary cosmological constant, T being the scale. 
If X > - VA then X(T) ^ v^A for T ^ oo, VA being a "fixed point" (a zero of 
the "/5-function" -(X^ - A) in eq. (4.3)). 

Let Li denote the length of the entry boundary and L2 the length of the exit 
boundary. Rather than consider a situation where the boundary cosmological 
constant X is fixed we can consider Li as fixed. We denote the corresponding 
propagator Ga{Li, Y; T). Similarly we can define G'a(X, L2] T) and Ga{Li, L2; T). 
They are related by Laplace transformations. For instance: 

Ga{X, Y-T)= / dL2 / dLi G(Li, L2; T) e-^^l-^^^ (4.5) 

JQ Jo 
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and one has the following composition rule for the propagator: 

POO 

GA(X,r;Ti + T2)= / dLGA{X,L;T,)G{L,Y,T2). (4.6) 

Jo 

We can now calculate the expectation value of the length of the spatial slice at 
proper time t G [0, T]: 

1 

W))x,Y,T= ^^^^ y.^^ j^ dLGA{X,L;t)LGA{L,Y;T-t). (4.7) 

In general there is no reason to expect {L{t)) to have a classical limit. Consider for 
instance the situation where X and Y are larger than and where T ^ 1/ -\/A. 
The average boundary lengths will be of order 1/X and 1/Y. But for ^ t ^ T 
the system has forgotten everything about the boundaries and the expectation 
value of L(t) is, up to corrections of order e~^^* or e-2^/A(T-^)^ determined by the 
ground state of the effective Hamiltonian Heff corresponding to Ga(X, Y; T). One 
finds for this ground state (L) = l/y/X. This picture is confirmed by an explicit 
calculation using eq. (4.7) as long as X, > V^. The system is thus, except for 
boundary effects, entirely determined by the quantum fluctuations of the ground 
state of Heff. 

We will here be interested in a different and more interesting situation where a 
non-compact spacetime is obtained as a limit of the compact spacetime described 
by (4.7). Thus we want to take T (yo and at the same time also the length 
of the boundary corresponding to proper time T to infinity. Since T — oo forces 
X(T, X) — >■ a/a it follows from (4.2) that the only choice of boundary cosmological 
constant Y independent of T, where the length {L{T)) x,y,t goes to infinity for 
T — >■ oo is y=— (fig. 4.1), since we have: 

(Tm\ - 1 gGA(X,y;T) _ 1 

With the choice Y = —\^ one obtains from (4.7) in the limit T ^ oo: 
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{Lit))x = ^ sinh (2v^(i + toiX))^ , (4.9) 
where to{X) is defined in eq. (4.4). 

We have caUed L2 the (spatial) length of the boundary corresponding to T and 
{L{t))x the spatial length of a time-slice at time t in order to be in accordance with 
earlier notation [30, 21], but starting from a lattice regularization and taking the 
continuum limit L is only determined up to a constant of proportionality which we 
fix by comparing with a continuum effective action. In section (2.7.3) wc showed 
that such a comparison leads to the identification of L as Lcont/T^ and we are led 
to the following 

L^t{t) = n{L{t))x = sinh (2^(t + to(X))) . (4.10) 



A 

Consider the classical surface where the intrinsic geometry is defined by proper 
time t and spatial length Lcontit) of the curve corresponding to constant t. It has 
the line element 

r2 smh^ (2VA{t + to{X))) 

ds' = + =^ de^ = dt^ + ^ ^ de^, (4.11) 

47r2 4A 

where t > and to{X) is a function of the boundary cosmological constant X at 
the boundary corresponding to i = (see eq. (4.4)). What is remarkable about 
formula (4.11) is that the surfaces for different boundary cosmological constants 
X can be viewed as part of the same surface, the Poincarc disc with curvature 
R = — 8A, since t can be continued to t = —tg. The Poincare disc itself is formally 
obtained in the limit X — > 00 since an infinite boundary cosmological constant 
will contract the boundary to a point. 



4.3 The classical effective action 

In this section we make a small digression to the "classical" theory and show that 
the emergence of the hyperbolic plane is natural from this point of view. In section 
(2.7.3) we discussed the derivation of the quantum Hamiltonian of causal quantum 
gravity from the following classical action. 

To make contact with the inherently quantum calculation by causal dynamical 
triangulations it is interesting to look at the classical behavior corresponding to 
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this action. The classical solutions corresponding to action (4.12) are 

l{t) = ^ sinh2VAt, k > elliptic case, (4.13) 
v A 

l{t) — "^C^- cosh2-\/Ai, K<0 hyperbolic case, (4-14) 
VA 

lit) = e^^*, K = parabolic case, (4.15) 

all corresponding to cylinders with constant negative curvature — 8A. In the elliptic 
case, where t must be larger than zero, there is a conical singularity at t = unless 
For n — \ the geometry is regular at t = and this value of k corresponds 
precisely to the Poincare disc, t — ^ being the "center" of the disc. So we see that 
for 1 the classical solution coincides nicely with the emergent geometry derived 
in the previous section 



4.4 Quantum fluctuations 

In many ways it is more natural to fix the boundary cosmological constant than to 
fix the length of the boundary. However, one pays the price that the fluctuations 
of the boundary size are large, in fact of the order of the average length of the 
boundary itself ^: from (4.8) we have 

(L^(r))x,y;T - {L{T))\y.^^ = jm))xy-,T ^ ^^(^r))^,^^^. (4.16) 

Such large fluctuations are also present around {L{t))xy;T for t < T. Prom this 
point of view it is even more remarkable that {L{t)) y=-^-t=oo ^^^^ ^ ^^^^ 
semiclassical interpretation. Let us now by hand fix the boundary lengths Li 
and L2. This is done in the Hartle-Hawking Euclidean path integral when the 
geometries [g] are fixed at the boundaries [68] . For our one-dimensional boundaries 
the geometries at the boundaries are uniquely fixed by specifying the lengths of 
the boundaries, and the relation between the propagator with fixed boundary 
cosmological constants and with fixed boundary lengths is given by a Laplace 
transformation as shown in eq. (4.5). Let us for simplicity analyze the situation 
where we take the length Li of the entrance loop to zero by taking the boundary 
cosmological constant X — > 00. Using the decomposition property (4.6) one can 
calculate the connected "loop- loop" correlator for fixed L2 and < t < t + A < T, 

{L{t)L{t + A))g,T = {L{t + A)L{t))L,,T - {L{t)){L{t + A))^„t. (4.17) 
One finds 

^This is true also in Liouville quantum theory, the derivation is essentially the same as that 
given in (4.16), as is clear from [66]. 
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2 sinh^ V^tsinh^ V^(T-(t+A)) 
A sinh^ 
2L2 sinh^ VAt sinh VA(t + A) sinh VA(T - (t + A) ) 

Va sinh^ Vat 



We also note that 



2 sinh VAi sinh VA(T-t) ^ sinh^ /At 

For fixed L2 and T — > 00 we obtain 

{L{t)L{t + A))g = e-^^^ (1 - e-^^*) ' (4.20) 



and 



m)L, = (1 - e-^v^*) . (4.21) 



Eqs. (4.20) and (4.21) tell us that except for small t we have {L{t))L2 = 1/\/A- 
The quantum fluctuations AL{t) of L{t) are deflned by {AL{t)y = (L(t)L(t))('^). 
Thus the spatial extension of the universe is just quantum size (i.e. l/\/A, A 
being the only coupling constant) with fluctuations AL{t) of the same size. The 
time correlation between L{t) and L{t + A) is also dictated by the scale 1/\/A, 
telling us that the correlation between spatial elements of size 1 / \/A, separated in 
time by A falls of exponentially as e"^^^ . The above picture is precisely what 
one would expect from the classical action, which is proportional to the area and 
the boundary cosmological constants only, if we force T to be large and choose a 
Y such that (L2(T)) is not large, the universe will be a thin tube, "classically" of 
zero width, but due to quantum fluctuations of average width 1/VX. 
A more interesting situation is obtained if we choose Y — —y/X, the special value 
needed to obtain a non-compact geometry in the limit T — > 00. To implement this 
in a setting where L2 is not allowed to fluctuate we flx L2{T) to the average value 
(4.8) for r=-\/A: 

^^m-i^m).y^-V-.r-^, ^,,J^^_, (4.22) 
Prom (4.18) and (4.19) we have in the hmit T — > 00: 

{L{t)) = sinh2VVVi, (4.23) 
V A 



in accordance with (4.9), and for the "loop- loop" -correlator 
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{L{t + A)L{t))(^^ = I sinh^ "^^^ix (^^^^^^ ~7xi^~ ■ ^^-2^^ 

It is seen that the "loop-loop" -correlator is independent of A. In particular we 
have for A = 0: 



for t ^ 1/a/A- The interpretation of eq. (4.25) is in accordance with the picture 
presented below (4.21): We can view the curve of length L{t) as consisting of 
N{t) ~ \/AL(t) pa e^^* independently fluctuating parts of size and each 

with a fluctuation of size I/VA- Thus the total fluctuation AL{t) of L{t) will be 
of order 1/VA x y/N{t), 



VX{L{t)) 



~ e-^^*, (4.26) 



i.e. the fluctuation of L{t) around {L{t)) is small for t » I/a/A- In the same 
way the independence of the "loop- loop" -correlator of A can be understood as the 
combined result of L(t + A) growing exponentially in length with a factor e^^'^ 
compared to L{t) and, according to (4.20), the correlation of "line-elements" of 
L{t) and L{t + A) decreasing by a factor e"^"^^. 



4.5 Summary 

We have described how the CDT quantization of 2d gravity for a special value of 
the boundary cosmological constant leads to a non-compact (Euclidean) AdS-like 

spacetime of constant negative curvature dressed with quantum fluctuations. It is 
possible to achieve this non-compact geometry as a limit of a compact geometry 
as described above. In particular the assignment (4.22) leads to a simple picture 
where the fluctuation of L{t) is small compared to the average value of L{t). In 
fact the geometry can be viewed as that of the Poincare disc with fluctuations 
correlated only over a distance 1 / V^- 

Our construction is similar to the analysis of ZZ-branes appearing as a limit of 
compact 2d geometries in Liouville quantum gravity [66]. In the CDT case the non- 
compactness came when the running boundary cosmological constant X{T) went 
to the flxed point y/X for T — > oo. In the case of Liouville gravity, represented by 

DT (or equivalently matrix models), the non-compactness arose when the running 
(Liouville) boundary cosmological constant XLiouvUieiT) went to the value where 

the disc amplitude PF^(X) = 0, i.e. to X = \/a/2 (see eq. (4.1)). It is the same 
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process in the two cases, since the relation between Liouville gravity and CDT is 
well established and summarized by the mapping [31]: 




(4.27) 



between the coupling constants of the two theories. The physical interpretation 
of this relation is discussed in [31, 21]: One obtains the CDT model by chopping 
away all baby-universes from the Liouville gravity theory, i.e. universes connected 
to the "parent-universe" by a worm-hole of cut-off scale, and this produces the 
relation (4.27) ^. It is seen that X \/K corresponds precisely to X — > \f\j2. 
While the starting point of the CDT quantization was the desire to include only 
Lorentzian, causal geometries in the path integral, the result (4.11) shows that 
after rotation to Euclidean signature this prescription is in a natural correspon- 
dence with the Euclidean Hartle- Hawking no-boundary condition, since all of the 
geometries (4.11) have a continuation to t = — to, where the spacetime is regular. It 
would be interesting if this could be promoted to a general principle also in higher 
dimensions. The computer simulations reported in [24, 25, 26, 27, 28, 29] seem in 
accordance with this possibility. 



^The relation (4.27) is similar to the one encountered in regularized bosonic string theory in 
dimensions d > 2 [69, 44, 70]: The world sheet degenerates into so-called branched polymer. 
The two-point function of these branched polymers is related to the ordinary two-point function 
of the free relativistic particle by chopping off (i.e. integrating out) the branches, just leaving 
for each branched polymer connecting two points in target space one jyath connecting the two 
points. The mass-parameter of the particle is then related to the corresponding parameter in the 

partition function for the branched polymers as Xjs/H to X/\/L. in eq. (4.27). 
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Topology fluctuations of space 

and time 



In this chapter we refocus our attention on the issue of topology change in quantum 
gravity. We make some steps to go beyond the results of chapter 3 and incorporate 
topology fluctuations of space and time in the path integral of two dimensional 
CDT. In field theory language one would say that in this chapter we go beyond 
tree level by computing loop corrections. The incorporation of loop corrections 
is a notoriously difficult undertaking however. Even in standard quantum field 
theory the loop expansion is not convergent so one can only validly utilize the 
perturbation theory up to a limited number of loops. Technically, the expansion 
merely forms an asymptotic power series instead of a convergent one. The basic 
reason behind the divergence of the series expansion is that for large number of 
loops the amount of diagrams grows super-exponentially. 

The situation for gravity is typically similar, the number of ways in which one 
can cut and reglue a manifold to obtain manifolds with a different topology is 
very large and leads to uncontrollable divergences in the path integral. Even in 
the simplest case of two dimensional Euclidean geometries, the number of possible 
configurations grows faster than exponentially with the volume of the geometry. 
A well-known manifestation of this problem is the non-Borel summability of the 
genus expansion in string theory. This does not necessarily mean that there is no 
underlying well-defined theory, but even in the much-studied case of two dimen- 
sional Euclidean quantum gravity no physically satisfactory, unambiguous solution 
has been found [18]. 

Our contribution to the problem of spacetime topology change is twofold. In 
section (5.1) we take a modest point of view and analyze the issue of spacetime 
topology change perturbatively in the loop expansion. Particularly, we expand 
the standard formalism of two dimensional causal quantum gravity and compute 
the Hartle Hawking wavefunction up to second order in the genus expansion. The 
statistical mechanics of loop diagrams is considerably more complicated than the 
tree diagrams that were analyzed in chapter 3 henceforh we are not be able to 
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address the summability of the expansion. 

On a more positive and perhaps more speculative note we introduce a toy model in 
section (5.2) where we show that it might be possible to define a nonperturbative 
sum over topologies provided suitable causality restrictions are imposed ^. This 
section is based on [72]. 



5.1 Perturbation theory 

Recall that the time dependent disc amplitude of pure CDT is given by 

Wa^T) = — -^^e-^^-**^^^ (5.1) 
smh vAT 

which is obtained from the propagator by shrinking the initial boundary to a 
point, i.e. Wa{L, T) = G\{Li = 0, L; T). Upon integrating over time one obtains 
the Hartle Hawking wavefunction, 

W^{L) = j^ dTWj,(L,T) = ^—. (5.2) 

To go beyond this genus zero result we need to extend the existing literature on two 
dimensional CDT by computing amplitudes where the spatial topology is allowed 
to change. A simple example of such an object is the before mentioned trouser 
amphtude (fig. 5.1). It can be obtained by "gluing" three cyhnder amplitudes. 
In this gluing procedure one takes a propagator with a final loop length (L + L') 
that is equal to the sum of the lengths of the initial loops L and L' of two other 
propagators and integrates over L and L'. As in the case of the composition 
rule (2.32) one has to include a measure factor to obtain the correct amplitudes. 
Furthermore, from a continuum perspective one can say that the gluing is defined 
such that Dehn twists around spatial slices are not present. Instead of writing out 
the measure factors explicitly we absorb them into the propagators by marking 
the loops of the amplitudes that are being glued. For convenience we also use a 
mixed representation for the propagators where one of the boundaries has fixed 
cosmological constant and the other has fixed boundary length, 

-Xx{T)L -LXoo(T) 

G^(X,L,T)^^-^ ^ (5.3) 

where Xx{T) is given by 

Xx(T) = Xoc(T) , ^ ^ (5.4) 

^ ' ^ ^ sinh2VAT(X + Xoo(T)) ^ ^ 

and (T) = y/X coth y/AT. 

^Also in the context of two dimensional Euclidean quantum gravity a nonperturbative sum 
over genera has been performed in a simplified model [71]. 
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L+L 



X 



Figure 5.1: Elementary trouser amplitude by gluing propagators 



We can now calculate the trouser amplitude from an initial boundary with bound- 
ary cosmological constant X to two final boundaries of length L\ and L2 in time 
T, where the splitting occurs after a fixed time t (fig. 5.1), 



TA{X,Li,L2;T,t) = 2j j dLdL'Gf''^\x, L + L' ,t) x 



Gi^^°\L, Li, T-t)4^^°\L', L2, T-t), (5.5) 



where the superscript notation of the propagators denotes the marking of its loops, 
as introduced in section (2.7.2). Performing the integrations yields 



TA{X,Li,L2;T,t) = 2 6'^''^'^'^^^'^^^^ 



A cosh VAt + X sinh 
VA cosh Vat + X sinh VAT 



+ 2 e"^°°(*)(^i+^2) 



sinh ^/lit 
sinh 



(5.6) 



In the rest of the section we only present results where the length of the initial loop 
is zero to keep the presentation as transparent as possible. The trouser amplitude 
with zero initial length is dubbed the splitting amplitude which we denote by S. In 
particular we are interested in the splitting amplitude to two final loops of length 
Li and L2, where the time before the splitting is arbitrary and the time after the 
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splitting is fixed to be t' 

oo oo 

S^{L„W,t') = 2 j dL j dL'w'l\L + L')G^^'\L,L,-t')G^^'\L\L,-,t') 



= 2e-(^^+^')^e-^*'^. (5.7) 



This amplitude might seem a bit unnatural, since it treats the time intervals before 
and after splitting differently. It is however useful for studying averages of the 
duration of a hole (section (5.1.1)). If we integrate over the time after the split in 
spatial topology as well we obtain 

oo , 

r ■~{L1+L2)VA 

5a(Li,L2) = J dt'SA{L^,L2;t') = ^ . (5.8) 



In the following we use the results obtained in this subsection to calculate higher 
genus Hartle Hawking wavefunctions within our model. 



5.1.1 Higher genus Hartle Hawking wavefunctions 

In this section we present genus one and genus two generalizations of the Hartle 
Hawking wave function where we integrate over the time intervals. It is possible 

to obtain results where the time intervals are fixed, but we choose to omit them 
for sake of clarity. Given the splitting amplitude (5.8) it is straightforward to 
compute the genus one generalization of the Hartle Hawking wavefunction (5.2). 
To simplify calculations we first obtain the genus one wave function with fixed 
boundary cosmological constant Y and a mark on the boundary (see fig. 5.2), 



OO oo 

__^^Y)=e-^^ J dL J dL'Si^''\L,L')G^^-^'\L + L',Y), (5.9) 



where K = 2n/G'}^ is proportional to the inverse bare Newton's constant and S^^'^^ 
denotes that both loops of the splitting amplitude possess a mark. Since 2D 
quantum gravity can be viewed as string theory with a zero dimensional target 
space, e^'^ can be identified with a bare string coupling Qg. The renormalization 
of this coupling is addressed in the next section. Performing the integrations in 
(5.9) one obtains 

<..,K) = ,2r!±l:^f!±ii^:±i^. (5.10) 



32(^F + \/a) A5/2 
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If we now do an inverse Laplace transformation from Y to L and divide by L to 
remove the mark, we obtain the desired genus one Hartle Hawking wave function, 

WA,,=iiL) = gl ^ L (5.11) 

One can use the previous results to compute simple observables. For example, 
using (5.7) we can obtain the average duration of the hole in the genus one disc 
amplitude: 



(^)hoie = J jdt't' jdlj dL'S'^''\L,L'-t')G^l''\L + L',Y) 

= (5.12) 

The size of the hole is determined by 1/ which is the only length scale in the 
model. Similarly we get for the fluctuations 



(At)hole = V <^')hole - (^)hole = (5-13) 

It is interesting to observe that (At)j^^[^ = (t)hoic- This reflects the fact that as for 
the spatial geometry, where we have (AL) ~ (L), the geometry of the holes is also 
purely governed by quantum fluctuations. 

The procedure to obtain the genus two wave function is analogous to the one out- 
lined above. In practice it involves doing much more tedious calculations however. 
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since one has to glue more and more propagators. This does not introduce any 
fundamental complications and does not give one much more physical insight so 
we only present the result of the computations. To obtain the genus two wave 
function one needs to add contributions from three different diagrams, 



(5.14) 



which evaluates to, 



'(3r+28Z6+150«^+570«'i+1575P+3150P+4725/+4725) 
92160A"/2 



(5.15) 



where the dimensionless variable Z = Lv A is the length of the boundary in units 

of vx. 



Double scaling limit 

Looking at the Hartle Hawking wave functions for genus g = 0, 1, 2 , i.e. Eqs. (5.2), 
(5.11) and (5.15), one can see that they all have different dimensions. Specifically, 
the g = 0, 1, 2 wave functions have dimension of A^/^, A~^/^ and A~^^/^ respectively. 
This would imply that the different genus contributions cannot be added to a 
single wave function. In fact one can add them if one takes into account their 
wave function renormalization factors that appear when taking the continuum 
limit of the discrete sums. The reason why those factors did not appear in the 
previous sections is that for the case of fixed topology one absorbs these factors in 
the boundary states of the theory. If one would have kept those factors from the 
outset the g = 0, 1, 2 wave functions would all be dimensionless and would behave 
as a\/A, (oa/A)"^ and (ca/A)"^^ respectively, where a is the cutoff of the theory 
induced by the discrete lattice spacing. 

To remove this cutoff dependence on the lattice spacing a one can do a renormal- 
ization of the string coupling such that all higher genus diagrams will contribute 
in the continuum limit. The simultaneous scaling of the string coupling and the 
cosmological constant is called the "double scaling" limit. In particular, we obtain 




Observe that this scaling limit is the same as the scaling we considered in chapter 
3. This implies that the spatial and spacetime topology changes are in fact part of 
one topological expansion. Using this scaling limit we can now derive the Hartle 
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Hawking wave function up to order two in the genus expansion, 




e-'[3f +281^^+1501^ +5701* +15751^ +31501'^ +4:7251+4725) 



+ 0{g%) 



92160 



(5.17) 



where we again used the dimensionless variable I — Ly/X. 



5.1.2 Summary 

We have described how to include manifolds of higher genus in the path integral 
of two dimensional CDT. One of the main results is the computation of the Hartle 
Hawking wave function up to two loops in the genus expansion. Generalization 
to higher genus is in principle straightforward, however, the calculations become 
more and more cumbersome. It is interesting to note that the Hartle Hawking 
wave functions in the framework of CDT are very similar to the Hartle Hawk- 
ing amplitudes in Euclidian Dynamical Triangulations [18]. Using the method of 
loop equations it is possible to also obtain higher genus results for 2D Euclidean 
quantum gravity [47, 48]. A detailed comparison of these results might lead to a 
better understanding of the relationship between Euclidean and Lorentzian quan- 
tum gravity. 
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5.2 Nonperturbative sum over topologies? 

As shown in previous work, there is a well-defined nonperturbative gravitational 
path integral including an explicit sum over topologies in the setting of causal 
dynamical triangulations in two dimensions. In this section we derive a complete 
analytical solution of the quantum continuum dynamics of this model, obtained 
uniquely by means of a double-scaling limit. We show that the presence of in- 
finitesimal wormholes leads to a decrease of the effective cosmological constant, 
reminiscent of the suppression mechanism considered by Coleman and others in 
the four-dimensional Euclidean path integral. Remarkably, in the continuum limit 
we obtain a finite spacetime density of microscopic wormholes without assum- 
ing fundamental discreteness. This shows that one can in principle make sense 
of a gravitational path integral which includes a sum over spacetime topologies, 
provided suitable causality restrictions are imposed on the path integral histories. 

5.2.1 Outline 

A new idea to tame the divergences associated with spacetime topology changes 
in the path integral was advanced in [73] and implemented in a model of two di- 
mensional nonperturbative Lorentzian quantum gravity. The idea is to include a 
sum over topologies, or over some subclass of topologies, in the state sum, but 
to restrict this class further by certain geometric (as opposed to topological) con- 
straints. These constraints involve the causal (and therefore Lorentzian) structure 
of the spacetimes and thus would have no analogue in a purely Euclidean formula- 
tion. In the concrete two dimensional model considered in [73], the path integral 
is taken over a geometrically distinguished class of spacetimes with arbitrary num- 
bers of "wormholes", which violate causality only relatively mildly (see also [74]). 
As a consequence, the nonperturbative path integral turns out to be well defined. 
This is an extension of the central idea of the approach of causal dynamical trian- 
gulations, namely, to use physically motivated causality restrictions to make the 
gravitational path integral better behaved (see [75] for a review) . 
In section (5.2.3), we will present a complete analytical solution of the statistical 
model of two dimensional Lorentzian random geometries introduced in [73], whose 
starting point is a regularized sum over causal triangulated geometries including 
a sum over topologies. For a given genus (i.e. number of (worm)holes in the 
spacetime) not all possible triangulated geometries are included in the sum, but 
only those which satisfy certain causality constraints. As shown in [73], this makes 
the statistical model well defined, and an unambiguous continuum limit is obtained 
by taking a suitable double-scaling limit of the two coupling constants of the model, 
the gravitational or Newton's constant and the cosmological constant. The double- 
scaling limit presented here differs from the one found in [73, 74], where only the 
partition function for a single spacetime strip was evaluated. We will show that 
when one includes the boundary lengths of the strip explicitly - as is necessary 
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to obtain the full spacctimc dynamics - the natural renormalization of Newton's 
constant involves the boundary "cosmological" coupling constants conjugate to the 
boundary lengths. Although the holes we include exist only for an infinitesimal 
time, and we do not keep track of them explicitly in the states of the Hilbert space, 
their integrated effect is manifest in the continuum Hamiltonian of the resulting 
gravity theory. As wc will see, their presence leads to an effective lowering of the 
cosmological constant and therefore represents a concrete and nonperturbative 
implementation of an idea much discussed in the late eighties in the context of 
the ill-defined continuum path integral formulation of Euclidean quantum gravity 
(see, for example, [76, 77]). 

The remainder of this section is structured as follows. In the next subsection, 
we briefly describe how a nonperturbative theory of two dimensional Lorentzian 
quantum gravity can be obtained by the method of causal dynamical triangulation 
(CDT), and how a sum over topologies can be included. For a more detailed 
account of the construction of topology-changing spacetimes and the geometric 
reasoning behind the causality constraints we refer the reader to [73, 74]. The 
main result of Sec. (5.2.3) is the computation of the Laplace transform of the 
one-step propagator of the discrete model for arbitrary boundary geometries. In 
Sec. (5.2.4) we make a scaling ansatz for the coupling constants and show that 
just one of the choices for the scaling of Newton's constant leads to a new and 
physically sensible continuum theory. We calculate the corresponding quantum 
Hamiltonian and its spectrum, as well as the full propagator of the theory. Using 
these results, wc compute several obscrvables of the continuum theory in Sec. 
(5.2.5), most importantly, the expectation value of the number of holes and its 
spacetime density. In Sec. (5.2.6), we summarize our results and draw a number 
of conclusions. In Appendix B, we discuss the properties of alternative scalings for 
Newton's constant which were discarded in the main text. This also establishes 
a connection with previous attempts [78, 79] to generalize the original Lorentzian 
model without topology changes. In Appendix C, we calculate the spacetime 
density of holes from a single infinitesimal spacetime strip. 

5.2.2 Implementing the sum over topologies 

Our aim is to calculate the (1+1) -dimensional gravitational path integral 



nonperturbatively by using the method of Causal Dynamical Triangulations (CDT). 
The sum in (5.18) denotes the inclusion in the path integral of a specific, causally 
preferred class of fiuctuations of the manifold topology. The action S{g^,y) consists 
of the usual Einstein-Hilbert curvature term and a cosmological constant term. 
Since we work in two dimensions, we recall that integrated curvature term is pro- 
portional to the Euler characteristic x = 2—20—6 of the spacetime manifold, where 




topol. 



(5.18) 
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Q denotes the genus (i.e. the number of handles or holes) and b the number of 
boundary components. Explicitly, the action reads 



where K — is the inverse Newton's constant and A the cosmological constant 

(with dimension of inverse length squared). 

Just like in the original CDT model [21] we will first regularize the path inte- 
gral (5.18) by a sum over piecewise flat two dimensional spacetimes, whose flat 
building blocks are identical Minkowskian triangles, all with one spacelike edge 
of squared length and two timclikc edges of squared length —aa^, where a 
is a real positive constant. The CDT path integral takes the form of a sum over 
triangulations, with each triangulation consisting of a sequence of spacetime strips 
of height At — 1 in the time direction. A single such strip is a set of h trian- 
gles pointing up and I2 triangles pointing down (fig. 2.7). Because the geometry 
has a sliced structure, one can easily Wick-rotate it to a triangulated manifold of 
Euclidean signature by analytically continuing the parameter a to a real negative 
value [80]. For simplicity, we will set a = —1 in evaluating the regularized, real 
and Wick-rotated version of the path integral (5.18). 

In the pure CDT model the one-dimensional spatial slices of constant proper time t 
are usually chosen as circles, resulting in cylindrical spacetime geometries. For our 
present purposes, we will enlarge this class of geometries by allowing the genus to be 
variable. We define the sum over topologies by performing surgery moves directly 
on the triangulations to obtain regularized versions of higher-genus manifolds [73, 
74]. They are generated by adding tiny wormholes that connect two regions of 
the same spacetime strip. Starting from a regular strip of topology [0, 1] x and 
height At = 1, one can construct a hole by identifying two of the strip's timelike 
edges and subsequently cutting open the geometry along this edge (fig. 5.3). By 
applying this procedure repeatedly (obeying certain causality constraints [73, 74]), 
more and more wormholes can be created. Once the regularized path integral 
has been performed, including a sum over geometries with wormholes, one takes 
a continuum limit by letting a — * and renormalizing the coupling constants 
appropriately, as will be described in the following sections. 

Note that our wormholes are minimally causality- and locality-violating in that 
they are located within a single proper-time step (the smallest time unit available 
in the discretized theory) and the associated baby universes which are born at time 
t are reglued at time t + 1 "without twist" [73, 74]. In a macroscopic interpretation 
one could describe them as wormholes which are instantaneous in the proper-time 
frame of an ensemble of freely falling observers. Note that this is invariantly 
defined (on Minkowski space, say) once an initial surface has been chosen. Such 
a restriction is necessary if one wants to arrive at a well-defined unitary evolution 
via a transfer matrix formalism, as we are doing. To include wormholes whose 
ends lie on different proper-time slices, one would have to invoke a third-quantized 




(5.19) 
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formulation, which would very likely result in macroscopic violations of causality, 
locality and therefore unitarity, something we are trying to avoid in the present 
model. 

One could wonder whether the effect in the continuum theory of our choice of 
wormholes is to single out a preferred coordinate system. Our final result will show 

that this is not the case, at least not over and above that of the pure gravitational 
model without topology changes. The effect of the inclusion of wormholes turns 
out to be a rather mild "dressing" of the original theory without holes. We believe 
that the essence of our model lies not so much in how the wormholes are connected, 
because they do not themselves acquire a nontrivial dynamics in the continuum 
limit. Rather, it is important that their number is sufficiently large to have an 
effect on the underlying geometry, but on the other hand sufficiently controlled so 
as to render the model computable. 

A similar type of wormhole has played a prominent role in past attempts to devise 
a mechanism to explain the smallness of the cosmological constant in the Euclidean 
path integral formulation of /oitr-dimensional quantum gravity in the continuum 
[76, 77]. The wormholes considered there resemble those of our toy model in that 
both arc non-local identifications of the spacctimc geometry of infinitesimal size. 
The counting of our wormholes is of course different since we are working in a 
genuinely Lorentzian setup where certain causality conditions have to be fulfilled. 
This enables us to do the sum over topologies completely explicitly. Whether a 
similar construction is possible also in higher dimensions is an interesting, but at 
this stage open question. 

5.2.3 Discrete solution: the one-step propagator 

For the (1+1) -dimensional Lorentzian gravity model including a sum over topolo- 
gies, the partition function of a single spacetime strip of infinitesimal duration 
with summed-over boundaries was evaluated in [73] and [74]. In the present sec- 
tion, we will extend this treatment by calculating the full one-step propagator, 
or, equivalently, the generating function for the partition function of a single strip 
with given, fixed boundary lengths. This opens the way for investigating the full 
dynamics of the model. 

The discrete set-up described above leads to the Wick-rotated one-step propagator 
GUh, k, i = 1) = e-^('^+'^) J2 (5-20) 

T\h,l2 

where k. is the bare inverse Newton's constant and A the bare (dimensionless) 
cosmological constant, and we have omitted an overall constant coming from the 
Gauss-Bonnet integration. The sum in (5.20) is to be taken over all triangulations 
with li spacelike links in the initial and I2 spacelike links in the final boundary. Note 
that the number of holes does not appear as one of the arguments of the one-step 
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Figure 5.3: Construction of a wormhole by identifying two timelike edges of a 
spacetime strip and cutting open the geometry along the edge. 

propagator since we only consider holes that exist within one strip. Consequently, 
the number of holes does not appear explicitly as label for the quantum states, 
and the Hilbert space coincides with that of the pure CDT model. Nevertheless, 
the integrated effect of the topologically non-trivial configurations changes the 
dynamics and the quantum Hamiltonian, as we shall sec. 
The one-step propagator (5.20) defines a transfer matrix T by 

GxAh,l2,l)^{l2\f\h), (5.21) 

from which we obtain the propagator for t time steps as usual by iteration, 

Gx,.{h,k,t)^{l2\f%). (5.22) 

For simplicity we perform the sum in (5.20) over triangulated strips with periodi- 
cally identified boundaries in the spatial direction and one marked timelike edge. 
By virtue of the latter, Ga,k(^i, ^2, ^) satisfies the desired composition property of 
a propagator, 

G'a,«;(^1, ^2, ^1 + ^2) — ^ ^^,^(^1, il)GA,K(^, ^2, ^2), (5.23) 

I 

Gx,.{hj2:t+1) = ^GA,«(/i,/,l)GA,«(i,i2,i), (5.24) 

I 

where the sums on the right-hand sides are performed over an intermediate constant- 
time slice of arbitrary discrete length I. 

Performing the fixed-genus part of the sum over triangulations in (5.20) yields 

«U...) = e-^"EO(fJ^e--. (...) 
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with N = I1 + I2. To simplify calculations we will use the generating function 
formalism with 

00 

G{x,y,g,h,l)= J2 GxAh,l2,l)xW\ (5-26) 

h,l2=0 

where we have defined g = and h = e~'^. The quantities x and y can be 
seen as purely technical devices, or alternatively as exponentiated bare boundary 
cosmological constants 

X = e~^'", y = e"^""*. (5.27) 

Upon evaluating the sum over li and I2 one obtains the generating function of the 
one-step propagator 

Gix, y, g, h, 1) = 1 ^ . , /^^ ' (5-28) 

1-g {x + y)l + Vl - 4x^2 

with 

/ ^ - (5.29) 

Note that in order to arrive at the final result (5.28), we have performed an explicit 
sum over all topologies! The fact that this infinite sum converges for appropriate 
values of the bare couplings has to do with the causality constraints imposed on the 
model, which were geometrically motivated in [73], and which effectively reduce 
the number of geometries in the genus expansion. 

In (5.28) one recognizes the generating function Catlu^) for the Catalan numbers, 

C^t{u') = L^ . (5.30) 
1 + Vl - 4-^2 

For h — one has Cat{u'^) = 1 and expression (5.28) reduces to the one-step 
propagator without topology changes, 

G{x, y,g,h = 0, 1) = ^— -. (5.31) 

1-g [x + y) 

Furthermore, one recovers the one-step partition function with summed-over bound- 
aries of [73, 74] by setting x — y — 1, 

Z{g, h, 1) = -\- ^ ■ (5.32) 

1-2^1 + ^1-4(3^^)2 

5.2.4 Taking the continuum limit 

Taking the continuum limit in the case without topology changes is fairly straight- 
forward [21]. The joint region of convergence of (5.31) is given by 

|a;|<l, \y\<l, \g\ < ^. (5.33) 
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One then tunes the couphngs to their critical values according to the scaling rela- 
tions 

g = ^{l-a'A) + Oia'), (5.34) 
x = l-aX + 0{a'^), y = l-aY + 0{a^). (5.35) 

Up to additive renormalizations, x, y and A scale canonically, with corresponding 
renormalized couplings X, Y and A. In the case with topology change we have to 
introduce an additional scaling relation for h. Since Newton's constant is dimen- 
sionless in two dimensions, there is no preferred canonical scaling for h. We make 
the multiplicative ansatz^ 

h ^ —=hren{adY, (5.36) 
v2 

where hren depends on the renormalized Newton's constant Gn according to 

hren ~ 6 ^/^'^ _ (5.37) 

In order to compensate the powers of the cut-off a in (5.36), d must have dimensions 
of inverse length. The most natural ansatz in terms of the dimensionful quantities 
available is 

d= (VA"(X + y)i-"). (5.38) 

The constants P and a in relations (5.36) and (5.38) must be chosen such as to 
obtain a physically sensible continuum theory. By this we mean that the one-step 
propagator should yield the Dirac delta-function to lowest order in a, and that 
the Hamiltonian should be bounded below and not depend on higher-order terms 
in (5.34), (5.35), in a way that would introduce a dependence on new couplings 
without an obvious physical interpretation. 

To calculate the Hamiltonian operator H we use the analogue of the composition 
law (5.23) for the Laplace transform of the one-step propagator [21], 

/dz 
-—G{x,z-';l)G{z,y,t). (5.39) 

In a similar manner we can write the time evolution of the wave function as 

When inserting the scaling relations (5.34), (5.35) and t—'^ into this equation it 
is convenient to treat separately the first factor in the one-step propagator (5.28), 
which is nothing but the one-step propagator without topology changes (5.31), 



^Here the factor is chosen to give a proper parametrization of the number of holes in terms 
of Newton's constant (see Section (5.2.5)). 



5.2 Nonperturbative sum over topologies? 



85 



and the second factor, the Catalan generating function (5.30). Expanding both 
sides of (5.40) to order a gives 

Oia^)) ,m . ^ { (^ + „ 1^) Ca.(«^) } HZ), 

(5.41) 

where we have used 

^(X,T + a) = e-"-^V(^,r), (5.42) 

with '0(X) = ■^(x = 1 — aX). Note that the first term on the right-hand side of 
(5.41), Laplace-transformed delta-function. The interesting new be- 

haviour of the Hamiltonian is contained in the expansion of the Catalan generating 
function. Combining (5.30) and (5.29), and inserting the scalings (5.34), (5.35), 
yields 

Cat(«2) = 1 + ^ren 2^3-2 ^ ^ (5 43^ 

{z-xf 

where h.o. refers to terms of higher order in a. In order to preserve the delta- 
function and have a non-vanishing contribution to the Hamiltonian one is thus 
naturally led to j3 = 3/2. For suitable choices of a it is also possible to obtain 
the delta-function by setting (3 = 1, but the resulting Hamiltonians turn out to be 
unphysical or at least do not have an interpretation as gravitational models with 
wormholes, as we will discuss in Appendix A.^ 
For P = 3/2 the right-hand side of (5.41) becomes 



dZ I 1 (2K-XZ 2 VA^" hi 



We observe that for a ^ the last term in (5.44) does not contribute to the 
Hamiltonian. Performing the integration for a > and discarding the possibility 
of fractional poles the Hamiltonian reads 

(5.45) 

For all a's, these Hamiltonians do not depend on higher-order terms in the scaling 
of the coupling constants. One can check this by explicitly introducing a term 
quadratic in a (which can potentially contribute to H) in the scaling relations 
(5.35), namely, 

X = l-aA: + ^70^X2 + 0(0^), 

y = 1-aY + ^ja'^Y^ + 0{a^), (5.46) 



^One might also consider scalings of the form h Cihren{0'd) + C2/iren(<it^)^^^i but they can 
be discarded by arguments similar to those of Appendix A. 
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and noticing that (5.45) does not depend on 7. After making an inverse Laplace 
transformation ip{L) = dXe^ ^iJj{X) to obtain a wave function in the "posi- 
tion" representation (where it depends on the spatial length L of the universe), 
and introducing m — 3a — 1 the Hamiltonian reads 

Since H is unbounded below for m ^ 2, we arc left with m = and m = 1 as 
possible choices for the scaling. However, setting m = merely has the effect of 
adding a constant term to the Hamiltonian, leading to a trivial phase factor for the 
wave function. We conclude that the only new and potentially interesting model 
corresponds to the scaling with m — l and 

h'^hl^A{X + Y)a', (5.48) 
with the Hamiltonian given by 

Note that for all values Gn > of the renormalized Newton's constant (5.36) the 
Hamiltonian is bounded from below and therefore well defined. It is self-adjoint 
with respect to the natural measure d^{L)—dL and has a discrete spectrum, with 
eigenfunctions 

^„(L) = Ae-V2A(i-'^?en)i'L,( 2v/2A(l-/i2^J), n = 0, 1, 2, ... , (5.50) 

where L„ denotes the n'th Laguerre polynomial. Choosing the normalization con- 
stants as 

An = ^8A(l-/i2^J, (5.51) 

the eigenvectors {^/^„(L), n = 0, 1,2, ...} form an orthonormal basis, and the corre- 
sponding eigenvalues are given by 

En = V2A(1 - hl,^) (2n + 1) , n = 0,l,2,.... (5.52) 

Having obtained the eigenvalues one can easily calculate the Euclidean partition 
function for finite time T (with time periodically identified) 



For completeness we also compute the finite-time propagator 

00 

Gk,gAL,,L2,T) = (L2|e-^^|Li) = ^ e-^^"V':(i^2)V'n(i^i). (5.54) 

n=0 
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Inserting (5.50) into (5.54) and using known relations for summing over Laguerre 
polynomials [81] yields 

smli(a;T) \ smli(a;r) / 

where we have used the shorthand notation uj — •y/2A(l — /i^en)- expected, for 
hren the results reduce to those of the pure two dimensional CDT model. 

5.2.5 Observables 

Due to the low dimensionality of our quantum-gravitational model, it has only a 
few observables which characterize its physical properties. Given the eigenfunc- 
tions (5.50) of the Hamiltonian (5.49) one can readily calculate the average spatial 
extension (L) of the universe and all higher moments 

/■oo 

(L-)„= / dLU'^\ML)\'. (5.56) 
Using integral relations for the Laguerre polynomials [81] one obtains^ 

m 

irmx f 1 Y V{n-m)V{m + l) 
\^ In - \ ^771 7^ ~F7:r"^vFr~z^ ^ 



^8A(l-/iLJ; r(n + l)r(-m) 
X 3F2(-n, 1 + m, 1 + m; 1, 1 + m - n; 1), (5.57) 
where 3^2(01, 02, 03; 61, 62; z) is the generalized hypergeometric function defined by 

3^2(01, 02, 03; 61, 62; ^) = > (u\ iu\ i\ ■ (5-58 

Observe that the moments scale as {L'^)^ ~ A~^ which indicates that the effective 
Hausdorff dimension is given by du = 2, just like in the pure CDT model [82]. 
In addition to these well-known geometric observables, the system possesses a new 
type of "topological" observable which involves the number of holes Ng, as already 
anticipated in [73, 74]. As spelled out there, the presence of holes in the quantum 
geometry and their density can be determined from light scattering. An interesting 
quantity to calculate is the average number of holes in a piece of spacetime of 
duration T, with initial and final spatial boundaries identified. Because of the 
simple dependence of the action on the genus this is easily computed by taking the 
derivative of the partition function Zt with respect to the corresponding coupling, 
namely, 

= ^^1^. (5.59) 



^Note that the poles of r(— m) cancel with those of the hypergeometric function. 
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Figure 5.4: The density of holes n in units of A as a function of Newton's constant 
Gn- 



Upon inserting (5.53) this yields 

coth (^2A{1 -h^.^jT 
In an analogous manner we can also calculate the average spacetime volume 

leading to 

{V) = T ^^^Z^-^ coth (vW^T^J t) . (5.62) 



^2A 

Dividing (5.60) by (5.62) we find that the spacetime density n of holes is constant, 

The density of holes in terms of the renormalized Newton's constant is given by 

n = ^ A. (5.64) 

gGjv _ 1 

The behaviour of n in terms of the renormalized Newton's constant is shown in 
fig. 5.4. The density of holes vanishes as Gn —>■ and the model reduces to the case 
without topology change. - An alternative calculation of the density of holes from 
an infinitesimal strip, which leads to the same result, is presented in Appendix B. 
We can now rewrite and interpret the Hamiltonian (5.49) in terms of physical 
quantities, namely, the cosmological scale A and the density of holes in units of A, 
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Figure 5.5: The coefficient of the effective potential, f{ri) = 1/(1 + ri), as function 
of the density of holes in units of A, rj = j-. 



i.e. ■'7 = -f, resulting in 



One sees explicitly that the topology fluctuations affect the dynamics since the 
effective potential depends on rj, as illustrated by fig. 5.5. 

It should be clear from expressions (5.65) and (5.64) that the model has two scales 
instead of the single one of the pure CDT model. As in the latter, the cosmological 
constant defines the global length scale of the two dimensional "universe" through 
(L) ~ The new scale in the model with topology change is the relative scale 
rj between the cosmological and topological fluctuations, which is parametrized by 
Newton's constant Gn- 

5.2.6 Summary 

In this section, we have presented the complete analytic solution of a previously 
proposed model [73] of two dimensional Lorentzian quantum gravity including a 
sum over topologies. The presence of causality constraints imposed on the path- 
integral histories - physically motivated in [73, 74] - enabled us to derive a new 
class of continuum theories by taking an unambiguously defined double-scaling 
limit of a statistical model of simplicially regularized spacetimes. After comput- 
ing the Laplace transform of the exact one-step propagator of the discrete model, 
we investigated a two-parameter family, defined by (5.36) and (5.38), of possible 
scalings for the gravitational (or Newton's) coupling, from which physical con- 
siderations singled out a unique one. For this case, we computed the quantum 
Hamiltonian, its spectrum and eigenfunctions, as well as the partition function 
and propagator. Using these continuum results, we then calculated a variety of 
physical observables, including the average spacetime density of holes and the 
expectation values of the spatial volume and all its moments. 




(5.65) 
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This should be contrasted with the previous treatment in [73, 74], in which only 
the one-step partition function with summed-over boundaries was evaluated. Be- 
cause of the lack of boundary information, no explicit Hamiltonian was obtained 
there. Moreover, it turns out that the scahng of the couphngs which in the cur- 
rent work led to the essentially unique Hamiltonian (5.49) could not have been 
obtained or even guessed in the previous work. This is simply a consequence of 
the fact that the dimensionful renormalized boundary cosmological constants make 
an explicit appearance in the scaling relation (5.48) for h, and thus for Newton's 
constant. We conclude that - unlike in the case of the original Lorentzian model 
- for two dimensional causal quantum gravity with topology changes one cannot 
obtain the correct scalings for the "bulk" coupling constants from the one-step 
partition function with boundaries summed over (which is easier to compute than 
the full one-step propagator). 

In contrast with what was extrapolated from the single-strip model in [74] , the to- 
tal number of holes in a finite patch of spacetime turns out to be a finite quantity 
determined by the cosmological and Newton's constants. Note that this finite- 
ness result has been obtained dynamically and without invoking any fundamental 
discreteness. Since the density of holes is finite and every hole in the model is 
infinitesimal, this implies - and is confirmed by explicit calculation - that the ex- 
pectation value of the number of holes in a general spatial slice of constant time 
is also infinitesimal. The fact that physically sensible observables are obtained in 
this toy model reiterates the earlier conclusion [73] that causality-inspired methods 
can be a useful tool in constructing gravitational path integrals which include a 
sum over topologies. 

Prom the effective potential displayed in fig. 5.5 one observes that the presence 
of wormholes in our model leads to a decrease of the "effective" cosmological 
constant f{ri)A. In Coleman's mechanism for driving the cosmological constant A 
to zero [76, 77], an additional sum over different baby universes is performed in 
the path integral, which leads to a distribution of the cosmological constant that is 
peaked near zero. We do not consider such an additional sum over baby universes, 
but instead have an explicit expression for the effective potential which shows 
that an increase in the number of wormholes is accompanied by a decrease of the 
"effective" cosmological constant. A first step in establishing whether an analogue 
of our suppression mechanism also exists in higher dimensions would be to try and 
understand whether one can identify a class of causally preferred topology changes 
which still leaves the sum over geometries exponentially bounded. 



6 

Conclusions 



Despite many attempts, gravity has resiliently resisted a unification with the laws 
of quantum mechanics. Besides a plethora of technical issues, one is also faced with 
many interesting conceptual problems. The study of quantum gravity in lower di- 
mensional models ameliorates the technical difficulties while still preserving some 
of the conceptually fascinating characteristics of quantum gravity. 

In this thesis we analyze the very simple model of two dimensional quantum grav- 
ity. Although a rather extreme simplification of four dimensional gravity, many of 
the most fundamental issues are still relevant. Moreover, two dimensional gravity 
is interesting since it can be viewed as a minimal version of string theory. 

The first fundamental aspect where we make a contribution is the problem of 
topology change of space. Particularly, we present an exactly solvable model which 
shows that it is possible to incorporate spatial topology changes in the path inte- 
gral rigorously. We show that if the change in spatial topology is accompanied by 
a coupling constant it is possible to evaluate the path integral to all orders in the 
coupling. Furthermore, the model can be viewed as a hybrid between causal and 
Euclidean dynamical triangulation models. An interesting avenue for further re- 
search is the question whether our model has an interpretation within string theory. 
In particular, is our new coupling constant really equivalent to the string coupling? 

The second conceptual topic that we cover is the emergence of geometry from 
a background independent path integral. We show that from a path integral 
over noncompact manifolds a classical geometry with constant negative curvature 
emerges. No initial singularity is present, so the model naturally realizes the Har- 
tle Hawking boundary condition. Furthermore, we demonstrate that under certain 
conditions the superimposed quantum fluctuations are small! The model is an 
interesting example where a classical background emerges from background inde- 
pendent quantum gravity in an exactly solvable setting. How does the emergent 
geometry behave? Can we make contact with effective descriptions of quantum 
geometry such as noncommutative geometry or doubly special relativity? The 
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answer to these questions is incomplete for now, but the exact solvabihty of the 
model suggests that at least a detailed analysis could be possible. 

To conclude, we tackle the problem of spacetime topology change. Although we 
are not able to completely solve the path integral over all manifolds with arbitrary 
topology, wc do obtain some results that indicate that such a path integral might 
be consistent if suitable causality restrictions are imposed. As a first step we ex- 
tend the existing formalism of causal dynamical triangulations by a perturbative 
computation of amplitudes that include manifolds up to genus two. Further we 
present a toy model where we make the approximation that the holes in the man- 
ifold are extremely small. This simplification allows us perform an explicit sum 
over all genera and analyze the continuum limit exactly. Remarkably, the pres- 
ence of the infinitesimal wormholes leads to a decrease in the effective cosmological 
constant, reminiscent of the suppression mechanism considered by Coleman and 
others in the four-dimensional Euclidean path integral. 

The results of this thesis show that we still know very little about the ultimate 
configmation space of quantum gravity. Even for the extremely simple case of two 
dimensional quantum gravity various new models can be constructed that seem to 
lead to a well defined theory of quantum geometry. To understand the situation 
even better it would be important to have a better continuum understanding of 
the results of causal dynamical triangulations. 



A 



Lorentzian triangles 



In this appendix, a brief summary of results on Lorentzian angles is presented, 
where we follow the treatment and conventions of [13]. 

Since in CDT one considers simplicial manifolds consisting of Minkowskian tri- 
angles, Lorentzian angles or "boosts" naturally appear in the Rcgge action as 
rotations around vertices. Recall from Section (2.5.3) that the definition of the 
scalar curvature at a vertex v is given by (2.25), 



where e„ = 27r — X^oi; ^« deficit angle at a vertex v and K is the dual volume 

of the vertex v. In general, the spacelike deficit angle can be positive or negative 
as illustrated in fig. A.l. Furthermore, if the deficit angle is timelike, as shown 
in fig. A.2, it will be complex. The timelike deficit angles are still additive, but 
contribute to the curvature (A.l) with the opposite sign. Hence, both spacelike 
defect and timelike excess increase the curvature, whereas spacelike excess and 
timelike defect decrease it. 

The complex nature of the timelikc deficit angles can be seen explicitly by noting 
that the angles 9i between two edges cij and 6j (as vectors in Minkowski space) are 
calculated using 



where (■,■) denotes the fiat Minkowskian scalar product and by definition, the 
square roots of negative arguments are positive imaginary. 

Having given a concrete meaning to Lorentzian angles, we can now use (A.2) to 
calculate the volume of Minkowskian triangles which we will then use to explicitly 
compute the volume terms of the Regge action. 

The triangulations we are considering consist of Minkowskian triangles with one 
spacelike edge of length squared = o? and two timelike edges of length squared 



— 277-, 



(A.l) 



cos Qi 




(A.2) 
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Figure A. 2: Illustration of a spacelike (a) and a timelike (b) Lorentzian deficit 
angle e^, at a vertex v. 



= —aa? with a > 0. The general argument a > is used to give a mathemat- 
ically precise prescription of the Wick rotation, but it can be set to a = 1 after 
the Wick rotation has been performed. With the use of (A. 2) we can calculate the 
volume of such a Minkowskian triangle, yielding 



Vol(triangle) = — V4a + 1. (A.3) 

Now one can define the Wick rotation W as the analytic continuation of a ^— —a 
through the lower-half plane. One then sees that for a > | under this prescription 
z Vol(triangle) i— >■ — Vol(triangle) (up to a 0{1) constant which can be absorbed in 
the corresponding coupling constant in the action). This ensures that 

In the following we set a = 1 again. Generalizations of this treatment to dimension 
(i = 3, 4 can be found in [80]. 



B 



Alternative scalings 



In this appendix we discuss the scahngs with [3 = 1 which we discarded as unphys- 
ical in Sec. (5.2.4) above. We proceed as before by inserting the scaling relations 
(5.34) and (5.46) into the composition law (5.40). Instead of using /3 = | we set 
/3 = 1, leading to the scahng 

h=^hrenaVr{X + Y)^-'', (B.l) 

where the normalization factor on the right-hand side has been chosen for later 
convenience. Up to first order in a one obtains 

/ioo 1^ 

— {A{X, Z) + B{X, Z)a + 0{a')} tp{Z), (B.2) 
-ioo '^'^^ 

where the leading-order contribution is given by 
with 

C{X, Z) = Vl - hl,^{X - Z)-2"A". (B.4) 

For the Laplace transform of A{X, Z) to yield a delta- function, the scaling should 
be chosen such that a ^ 0. Considering now the terms of first order in a, 



B{X,Z) = 



hlSX + Z - AZ^)k^ 



{X - zy+^"C{x, z) (1 + c(x, z)f 



XZ-2K + ,iX-Zf 



{x-zYC{x, z)(i + c(x, z)y 



one finds that for a ^ — 1 the continuum limit is independent of any "hole contri- 
bution" (i.e. terms depending on /i^en) and therefore leads to the usual Lorentzian 
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model. This becomes clear when one expands the last term of (B.5) in {X — Z), 
resulting in 

(B.6) 

For q; ^ — 1 the term depending on hren does not have a pole and therefore does 
not contribute to the Hamiltonian. Since we are only interested in non-fractional 
poles, this leaves as possible a- values only a — and a — 



B.l The case p=l, a = 

For a = the Hamiltonian retains a 7-dependence contained in the first line of 
(B.5). Since there is no immediate physical interpretation of 7 in our model, it 
seems natural to choose 7 = 0, although strictly speaking this does not resolve 
the problem of explaining the 7-dependence of the continuum limit. Setting this 
question aside, one may simply look at the resulting model as an interesting inte- 
grable model in its own right. In order to obtain a delta- function to leading order, 
one still needs to normalize the transfer matrix by a constant factor 2/(1 -|- s), 
with s :— a/I — /i^en- After setting 7 = and performing an inverse Laplace 
transformation, the Hamiltonian reads 

It is self-adjoint with respect to the measure dfi{L) = U~^dL. Further setting 

2 

L — ^ one encounters the one-dimensional Calogero Hamiltonian 

with u = and A = 1 — 4(1 — s)^, which implies that the model covers the 
parameter range — 3 ^ A ^ 1. The maximal range for which the Calogero Hamil- 
tonian is self-adjoint is —00 < A ^ 1. The usual Lorentzian model without holes 
corresponds to A — 1. The Hamiltonian ((B.8)) has already appeared in a causal 
dynamically triangulated model where the two dimensional geometries were dec- 
orated with a certain type of "outgrowth" or small "baby universes" [78]. This 
model covered the parameter range ^ A ^ 1. 
The eigenvectors of the Hamiltonian ((B.7)) are given by 

^n{L) = Ane""^^ iFi{-n, s, 2V2AL), (i/x(L) = L'-^dL, (B.9) 

where iFi{—n,a,b) is the Kummer confluent hypergeometric function. The eigen- 
vectors form an orthonormal basis with the normalization factors 
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and the corresponding eigenvalues 



v2A 

E^^- (2n + s), 71 = 0,1,2,.... (B.ll) 

s 

One sees explicitly that the case s = 1 or, equivalently, A = l corresponds to the 
pure two dimensional CDT model. 



B.2 The case = 1, a = 



For a — the result does not depend on 7 and therefore on the detailed manner 
in which we approach the critical point. However, the Hamiltonian 

cannot be made self-adjoint with respect to any measure dii{L) because the bound- 
ary part of the partial integration always gives a nonvanishing contribution. We 
therefore discard this possibihty. 



B Alternative scalings 



c 



The density of holes of an 
infinitesimal strip 



In this appendix we give an alternative derivation of the spacetime density n of 
holes and explicitly show that the number of holes in a spacetime strip of infinites- 
imal time duration a is also infinitesimal. The operator in the L-rcpresentation of 
the number of holes per infinitesimal strip with fixed initial boundary L can be 
calculated by 

N,,a-.0-f-^^^, (C.l) 

where T is the transfer matrix defined in (5.21). Using T — 1 — aH + O(a^) and 
evaluating (C.l) to leading order in a gives 

N,,a-.o --a^^ + Oia') = 2 A /^^ L a + 0{a% (C.2) 

Similarly, the volume operator of the same infinitesimal spacetime strip in the 
L-representation is given by 

K_o = -T-^— ^a— + 0{a') = 2 (1 - hlJLa + 0{a'). (C.3) 

Although both expressions (C.2) and (C.3) vanish in the limit as a — > (and 
therefore the number of holes and the strip volume are both "infinitesimal" ) , their 
quotient evaluates to a finite number independent of L, namely, 

^^iVg^^ n ^ ^^^^^ 

Va^O 1 — Alren 

This is the exactly the same result for the spacetime density n of holes as we 
obtained earlier from the continuum partition function (5.63). 
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